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SUMMARY 


Existing  numerical  models  of  river  channel  morphology  are  seriously  limited  in 
their  application  by  their  inability  to  account  for  bank  erosion  and  changing  channel  width. 
Recent  research  has,  however,  now  led  to  the  development  of  a  physically-based  numerical 
model  of  river  channel  widening,  in  which  changing  channel  width  is  accounted  for  by 
coupling  fluvial  bank  erosion  and  mass-wasting  algorithms  with  flow  and  sediment 
transport  models  (Darby,  1994;  Darby  and  Thome,  In  Review).  Verification  studies 
indicate  that  this  new  model  is  capable  of  predicting  channel  adjustments  well,  qualitatively, 
but  widening  rates  are  seriously  under-predicted  by  the  new  model  (Darby  et  ah.  In 
Review).  This  deficiency  has  focussed  attention  on  the  predictive  ability  of  versions  of  the 
Osman  and  Thome  (1988)  mass-wasting  algorithm  further  developed  by  Darby  (1994)  and 
currently  employed  in  the  Darby-Thorne  numerical  model.  In  particular,  the  Osman- 
Thome-Darby  analysis  constrains  the  failure  plane  to  pass  through  the  toe  of  the  bank, 
excluding  the  possibility  of  simulating  secondary,  upper  bank  failures.  The  Osman-Thome- 
Darby  analysis  also  excludes  pore  water  pressure  and  hydrostatic  confining  pressure 
effects.  These  deficiencies  have  been  addressed  conceptually  by  Simon  et  al.  (1991),  but 
this  approach  was  developed  only  for  banks  of  simple  geometry,  and  does  not  account  for 
the  effects  of  pore  and  confining  pressures  on  the  development  of  the  incipient  failure 
plane.  In  contrast,  the  Osman-Thome-Darby  analysis  utilises  a  more  complete,  physically- 
based  approach  in  predicting  failure  plane  angle,  and  the  approach  is  also  based  on  a  bank 
profile  characteristic  of  natural  eroding  riverbanks. 

In  this  report,  a  new  bank  stability  analysis  is  developed  which  combines  the 
advantages  of  the  Osman-Thome-Darby  analysis  with  the  advantages  of  the  Simon  et  al. 
(1991)  analysis,  with  the  aim  being  to  improve  the  predictive  ability  of  the  mass- wasting 
algorithm  and,  therefore,  the  predictive  ability  of  the  numerical  channel  widening  model. 
The  new  analysis  is  applicable  to  steep,  cohesive,  non-layered  riverbanks  which  fail  along 
planar  failure  surfaces.  The  analysis  is  developed  for  banks  with  profiles  characteristic  of 
natural,  eroding  riverbanks.  Failure  plane  angles  are  calculated  using  a  physically-based 
approach;  additionally  the  failure  plane  is  not  constrained  to  pass  through  the  toe  of  the 
bank.  Pore  water  and  hydrostatic  confining  pressure  effects  are  included  in  the  analysis. 
Using  the  procedure  developed  by  Darby  (1994),  the  new  analysis  also  is  capable  of 
predicting  the  probability  of  mass  failure,  as  well  as  the  overall  factor  of  safety  with  respect 
to  mass  failure. 

In  order  to  test  the  predictive  ability  of  the  new  analysis,  a  data  set,  comprising  51 
individual  mass-failures,  is  also  herein  compiled.  The  data  set  itself  is  in  a  spreadsheet 
format  on  a  diskette  supplied  with  this  report.  The  new  data  set  represents  the  first 
compilation  of  data  suitable  for  testing  the  ability  of  bank  stability  analyses  to  predict 
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accurately  the  factor  of  safety  and  geometry  of  individual  observed  mass-failures.  The  data 
set  has  been  used  to  compare  the  predictive  abilities  of  4  bank  stability  analyses:  Lohnes 
and  Handy  (1968),  Huang  (1983),  Osman  and  Thome  (1988)  and  the  analysis  developed 
in  this  report.  Results  indicate  that  the  analysis  developed  in  this  report  is  the  most 
successful  of  the  tested  analyses  in  predicting  the  stability  of  riverbanks  with  respect  to 
mass  failure.  However,  the  new  analysis  still  significantly  over-predicts  bank  stability. 
Nevertheless,  it  is  demonstrated  that  the  new  analysis  provides  significant  improvements  to 
predictions  of  bank  stability  when  compared  to  the  Osman  and  Thome  (1988)  analysis.  It 
is,  therefore,  likely  that  application  of  the  new  analysis  developed  herein  in  the  Darby- 
Thome  numerical  channel  evolution  model  would  result  in  improved  quantitative 
predictions  of  channel  adjustment.  Implementation  of  the  new  bank  stability  analysis  in  the 
numerical  channel  morphology  model  is  a  topic  of  current  research. 

In  the  meantime,  the  new  bank  stability  analysis  may  be  used  as  a  stand-alone  tool 
to  aid  engineers  in  the  design  of  stable  channels.  To  this  end,  the  new  analysis  has  been 
developed  as  a  computer  program,  which  is  supplied  on  disk  with  this  report,  together  with 
supporting  documentation.  The  computer  program  may  be  used  to  calculate  the  factor  of 
safety  or  probability  of  failure  of  riverbanks  of  geometry  specified  by  the  user.  The  user 
must  specify  the  geotechnical  characteristics  of  the  bank  material  and  select  whether  or  not 
to  include  pore  water  and/or  confining  pressure  effects  in  the  analysis.  If  the  user  does 
select  either  these  options  he  or  she  must  specify  the  ground  water  surface  elevation  and 
elevation  of  the  surface  water  in  the  channel,  respectively. 

A  particular  feature  of  the  computer  program  is  that  it  may  be  used  to  determine  the 
amounts  of  fluvial,  particle-by-particle,  erosion  and  bed  degradation  required  to  destabilize 
an  initially  stable  bank  to  the  point  of  failure.  This  feature  will  allow  the  computer  program 
to  be  utilised  as  a  design  aid  in  a  variety  of  stable  channel  design  applications,  such  as 
determining  minimum  grade-control  stracture  elevations  required  to  prevent  degradation 
initiating  channel  widening  and  instability  or  in  estimating  volumes  of  sediment  supplied  to 
the  channel  from  fluvial  erosion/bed  degradation  induced  bank  failures.  The  ability  of  the 
model  to  generate  predictions  of  probability  of  failure  will  also  allow  the  program  to  be 
used  as  a  tool  in  reliability  analyses  (e.g.  Johnson,  1992  ) 
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1.  INTRODUCTION 


Riverbank  erosion  and  associated  sedimentation  and  land  loss  hazards  are  a 
resource  management  problem  of  global  significance.  In  the  United  States,  an  estimated 
142,000  miles  of  stream  bank  are  in  need  of  erosion  protection,  with  the  cost  of  protecting 
US  stream  banks  in  1981  about  US$1  billion  annually  (U.S.  Array  Corps  of  Engineers, 
1983).  Not  only  may  infrastructure  located  near  destabilized  river  banks  be  threatened  by 
accelerated  bank  erosion  and  land  loss,  but  bank  erosion  often  generates  large  quantities  of 
fine-grained  sediments  that  may  pose  a  hazard  to  aquatic  habitats  as  well  as  contributing  to 
sedimentation  problems  downstream.  For  example,  the  U.  S.  Army  Corps  of  Engineers 
(1983)  study  of  the  sediment  budget  in  the  Sacramento  River  showed  that  of  the  11.5 
million  tons  of  total  sediment  inflow  to  the  system,  6.8  million  tons  (59%)  is  derived  from 
bank  erosion.  Such  river  bank  instability  is  often  triggered  by  human  modifications,  as  well 
as  by  longer  term  natural  modifications  caused  by  climate  change. 

Recently,  attempts  to  predict  the  land  loss  and  bank  sediment  yield  associated  with 
riverbank  instability  and  width  adjustments  (e.g.  Lohnes,  1991;  Darby  and  Thorne,  1992; 
1994)  have  been  made.  These  approaches  have  been  based  on  estimating  the  failure 
geometry  of  banks  that  become  unstable  following  bed  degradation  and/or  direct  fluvial 
erosion.  However,  the  predictive  ability  and  effectiveness  of  these  approaches  is  limited  by 
a  number  of  unrealistic  assumptions  used  in  the  analyses.  First,  the  geometry  of  the 
eroding  bank  is  usually  specified  in  an  unrealistic  way.  Most  stability  analyses  specify  the 
geometry  of  the  bank  profile  using  a  single  bank  height  and  angle.  In  fact  the  geometry  of 
natural  eroding  banks  is  more  complex,  due  to  processes  of  bed  degradation  or  aggradation 
and  direct  lateral  erosion  acting  to  deform  the  bank  face  (Osman  and  Thome,  1988).  Most 
simple  bank  stability  analyses  also  make  unnecessary  assumptions  constraining  the  failure 
surface  to  pass  through  the  toe  of  the  bank.  In  fact,  the  failure  surface  may  frequently  not 
pass  through  the  toe  (Simon  et  al.,  1991).  Finally,  the  influence  of  the  surface  and  ground- 
water  distributions  on  hydrostatic  confining  and  pore-water  pressures  acting  on  the 
incipient  failure  block  are  also  usually  neglected  (Simon  et  al.,  1991).  There  is  a  need  to 
address  these  limitations  in  order  to  improve  the  predictive  ability  of  methods  of  predicting 
land  loss  and  bank  sediment  yield. 

There  is  also  a  clear  need  to  predict  the  response  of  rivers  to  changes  in  control 
variables  and  river  engineering  works  in  order  to  assess  the  influence  of  environmental  and 
human-induced  changes  on  channel  form  and  also  because  channel  form  adjustments  and 
lateral  channel  migration  impact  flood  plain  dwellers  and  users  directly  through  erosion  and 
sedimentation. 

Studies  of  channel  morphology  have  demonstrated  the  significance  of  the  channel 
width  as  a  dominant  variable  in  controlling  process-form  interactions,  probably  through  the 
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influence  of  width  on  secondary  flow  scaling  (Yalin,  1971).  This  has  resulted  in  the 
widespread  use  of  the  channel  width  as  an  independent  or  specified  variable  in  the  various 
channel  morphology  models,  including  numerical  mobile-bed  models  such  as  HEC-6. 
HEC-6  is  the  preferred  one-dimensional  numerical  model  for  water  and  sediment  routing 
employed  by  Corps  of  Engineers  personnel  faced  with  real  world  problems.  However, 
approaches  which  use  the  channel  width  as  a  constant,  specified  variable  are  in  fact 
severely  limited  because  channel  changes  in  response  to  variations  in  flow  or  sediment 
regime  are  frequently  dominated  by  transient  width  adjustments  (Simon,  1992;  Darby, 
1994).  Consequently,  the  application  of  existing  numerical  models  of  channel  morphology, 
including  HEC-6,  is  limited  by  their  inability  to  account  for  bank  erosion  and  changing 
channel  width  through  time  (Darby,  1994).  The  neglect  of  a  widening  predictor  probably 
represents  the  most  significant  remaining  limitation  of  HEC-6,  and  inclusion  of  such  a 
predictor  would  represent  a  significant  advance  in  the  ability  of  HEC-6  to  model  stream 
channel  erosion  and  sedimentation  processes. 


Figure  1  Sediment  fluxes  in  the  near  bank  zone  (after  Thome  and  Osman.  19RR1 
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Recently,  advances  in  our  ability  to  understand  the  processes  and  mechanisms  of 
width  adjustment  in  river  channels  have  been  made  in  research  at  the  University  of 
Nottingham,  UK  (Darby,  1994;  Darby  and  Thome,  In  Review;  Darby  et  al..  In  Review). 
This  work  has  led  to  the  development  of  a  numerical  model  of  channel  moiphology  which 
includes  adjustment  of  the  channel  width  by  coupling  a  quasi  one-dimensional 
degradation/aggradation  model  with  a  bank  stability  algorithm  in  order  to  account  for 
temporal  trends  in  bank  geometry  and  stability  and  the  rate  of  retreat  of  the  banks.  This 
work  has  highlighted  the  importance  of  accurately  modelling  the  fluxes  of  sediment 
entering  and  leaving  the  near  bank  zone  (Figure  1)  in  order  to  model  the  rate  of  retreat  or 
advance  of  the  bank. 

The  research  at  Nottingham  has  shown  that  when  riverbanks  become  unstable  with 
respect  to  mass  failure  under  gravity,  it  is  the  inflow  of  bank  material  that  dominates  the 
near  bank  sediment  budget.  Consequently,  approaches  based  on  predicting  width  by 
coupling  bank  stability  models  with  bed  morphology  models  must  be  able  to  predict 
accurately  the  onset  of  riverbank  instability  and  the  magnitude  of  the  bank  sediment  flux, 
which  is  determined  by  the  geometry  of  the  mass  failure.  Indeed,  under-predictions  of 
widening  rates  in  verification  analyses  of  the  Darby  (1994)  model  were  attributed  primarily 
to  the  tendency  of  the  existing  (the  Osman-Thome-Darby)  bank  stability  algorithm  to 
under-predict  failure  block  widths  (Darby  et  al..  In  Review).  Development  of  a  reliable 
bank  stability  algorithm  is,  therefore,  also  fundamental  to  the  success  of  a  widening 
predictor  in  a  morphological  model  such  as  HEC-6.  The  development  of  an  improved  bank 
stability  algorithm  is  also  justified,  therefore,  through  the  potential  improvement  in  the 
performance  of  HEC-6  that  may  be  attained  by  incorporating  the  bank  stability  algorithm 
into  the  HEC-6  code. 

One  problem  in  evaluating  objectively  the  predictive  ability  of  riverbank  stability 
analyses  lies  in  the  difficulties  associated  with  measuring  the  geometries  of  eroding, 
unstable  banks  in  the  field.  This  has  meant  that  there  has  up  until  now  been  a  limited  data 
set  available  to  characterize  natural  unstable  riverbanks  at  the  point  of  failure.  This  has 
resulted  in  an  inability  to  assess  definitively  the  predictive  ability  of  the  wide  range  of 
existing  bank  stability  analyses,  leading  to  confusion  and  arbitrary  selection  of  analyses  by 
engineers  faced  with  real  world  problems.  It  is  apparent  that  collation  of  an  improved 
database  would  aid  greatly  predictive  comparisons  of  existing  and  new  analyses  of 
riverbank  stability,  helping  to  decide  which  stability  analyses  are  most  appropriate  for 
analysing  land  loss,  sediment  yield  impacts  and  channel  adjustments  resulting  from  bank 
erosion  processes. 


9 


2.  OBJECTIVES 


To  address  directly  the  problems  outlined  in  the  preceding  section,  the  objectives  of 
this  study  are  to: 

1.  Develop  an  improved  analysis  of  the  stability  of  natural,  eroding,  steep  cohesive 
riverbanks  which  fail  along  planar  surfaces; 

2.  Collate  a  set  of  data  suitable  for  testing  the  new,  and  existing,  analyses  of  bank  stability; 

3.  Determine  the  predictive  abilities  of  the  new  and  existing  analyses  of  riverbank  stability 
using  the  data  set  obtained  in  objective  (2); 

4.  Develop  a  computer  model  of  bank  stability  based  on  the  new  bank  stability  analysis. 
This  model  will  be  developed  bearing  in  mind  the  possibility  of  later  incorporating  the  new 
analysis  into  the  HEC-6  morphological  model. 

3.  RIVERBANK  STABILITY  ANALYSIS 

Over  the  years,  a  large  number  of  slope  stability  models  have  been  developed  which 
have  been  applied  to  analyse  the  stability  of  eroding  cohesive  riverbanks.  Reviews  of  these 
analyses  have  been  made  by  Thome  (1982)  and  Huang  (1983). 

Traditionally,  analyses  of  the  stability  of  steep,  cohesive  eroding  banks  which  fail 
along  planar  surfaces  have  been  based  on  estimating  the  driving  and  resisting  forces  acting 
on  incipient  failure  blocks  with  a  relatively  simple  geometry  (Figure  2),  in  which  the  failure 
block  is  described  in  terms  of  the  overall  bank  height  and  a  single  bank  angle  (that  is,  the 
bank  profile  has  only  one  surface).  Examples  of  analyses  based  on  geometries  of  this  type 
include  Lohnes  and  Handy  (1968),  Thorne  et  al  (1981)  and  Huang  (1983).  In  these 
analyses,  stability  is  modelled  by  defining  a  factor  of  safety  (FS)  as  the  ratio  of  resisting 
and  driving  forces  acting  on  the  incipient  failure  bock.  Failure  is  predicted  to  occur  when 
FS  <  I,  the  dimensions  of  the  failure  being  determined  from  the  dimensions  of  the  critical 
failure  surface  (Figure  2). 

However,  in  practice,  these  previous  analyses  are  limited  by  four  major 
shortcomings  when  applied  to  natural,  eroding  riverbanks.  First,  the  simple  geometry  used 
to  characterize  the  dimensions  of  the  incipient  failure  block  shown  in  Figure  2  is  inadequate 
to  describe  the  geometry  of  natural  eroding  riverbanks  (Osman  and  Thome,  1988).  This  is 
because  the  profile  of  natural  eroding  riverbanks  is  deformed  by  combinations  of  bed 
degradation  or  aggradation,  which  influence  the  overall  height  of  the  bank,  as  well  as  direct 
lateral  fluvial  erosion,  which  acts  to  steepen  the  bank  profile.  These  combinations  of  fluvial 
processes  tend  to  result  in  a  characteristic  geometry  of  natural  eroding  riverbanks  which 
may  be  more  realistically  illustrated  by  the  schematic  diagram  shown  in  Figure  3.  Second, 


10 


the  failure  plane  is  constrained  to  pass  through  the  toe  of  the  bank.  This  is  unrealistic,  since 
the  most  critical  failure  surface  may  intersect  the  bank  profile  at  any  point  on  the  bank 
profile  (Simon  et  al.,  1991).  A  third  limitation  is  that  the  influences  of  soil  pore  water 
pressures  and  the  hydrostatic  confining  pressure  of  surface  water  in  the  channel  adjacent  to 
the  bank  on  the  force  distribution  on  the  incipient  failure  block  is  either  usually  ignored,  or 
characterized  by  an  overly  simplistic  pore  pressure  ratio  term  (Simon  et  al,  1991).  Finally, 
the  failure  plane  angle  is  not  usually  estimated  using  rigorous  physically-based  methods, 
but  is  instead  frequently  approximated  as  a  function  of  bank  angle  and  soil  friction  angle 
only. 


A  B 


Many  of  these  limitations  have  individually,  or  in  combination,  been  addressed  in 
some  more  sophisticated  analyses.  Osman  and  Thome  (1988)  derived  an  analysis  of  bank 
stability  based  on  the  more  realistic  geometry  shown  in  Figure  3.  They  also  calculated 
failure  plane  angle  using  physically-based  techniques.  However,  pore  water  pressures  and 
hydrostatic  confining  pressures  were  neglected  in  their  analysis,  and  the  failure  plane  was 
constrained  to  pass  through  the  toe  of  the  bank.  Simon  et  al.  (1991)  developed  a  model  of 
bank  stability  which  did  account  for  hydrostatic  confining  and  pore  water  pressures,  and  in 


which  the  failure  plane  was  not  constrained  to  pass  through  the  toe.  However,  their  model 
was  based  on  the  simple  bank  geometry  of  Figure  2,  and  failure  plane  angles  were 
estimated  using  an  overly-simplistic  technique.  Furthermore,  there  are  some  deficiencies  in 
the  Simon  et  al.  analysis.  Notably,  this  analysis  does  not  account  for  the  influence  of 
dynamically  interacting  pore  and  surface  waters  on  the  pore  water  and  hydrostatic  pressure 
distributions.  No  model  has  yet  been  developed  which  attempts  to  address  each  of  the 
limitations  identified  above  in  a  single,  unified  analysis.  The  limitations  associated  with 
previous  analyses  should  mean  that  the  predictive  ability  of  these  approaches  is  not, 
therefore,  as  high  as  it  could  or  should  be. 


Lateral 

erosion 

Figure  3  Characteristic  geometry  of  natural,  eroding  riverbank 

In  response  to  the  need  to  predict  accurately  the  stability  and  geometry  of  incipient 
failure  blocks  along  unstable  river  channels,  a  new  analysis  is  developed  herein  in  which 
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each  of  these  limitations  are  addressed  specifically.  The  development  of  the  new  model  is 
now  described. 

3.1  Stability  Analysis 

The  basic  framework  for  the  bank  stability  analysis  is  illustrated  in  Figure  4,  which 
shows  a  bank  with  the  more  realistic  geometry  highlighted  by  Osman  and  Thome  (1988), 
but  in  which  the  failure  block  is  subjected  to  pore  water  pressures,  hydrostatic  confining 
pressures  and  in  which  the  failure  plane  is  free  to  intersect  the  bank  surface  at  any  point 
along  its  profile.  The  task  is  to  quantify  each  of  the  forces  acting  on  the  incipient  failure 
block,  including  the  effects  of  the  pore  water  and  hydrostatic  confining  pressures,  while 
still  allowing  for  the  failure  plane  to  intersect  the  bank  profile  at  the  most  critical  location. 
Additionally,  the  failure  plane  angle  should  be  estimated  using  physically-based  techniques. 


Figure  4  Bank  stability  analysis  definition  diagram 
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In  order  to  obtain  a  solution,  certain  assumptions  are  required  at  the  outset  of  the 
model  derivation.  These  assumptions  are  that  the  bank  is  non-layered  and  that  the  effects  of 
vegetation  on  bank  stability  may  be  excluded.  It  is  assumed  that  the  stability  of  the  bank 
may  be  modelled  using  the  factor  of  safety  concept,  where  failures  are  detected  once  the 
ratio  of  the  resisting  and  driving  forces  falls  below  unity.  The  factor  of  safety  is  defined  by: 


FS  = 


FR 

FD 


(1) 


where  FS  =  factor  of  safety,  FR  =  resultant  resisting  force  on  failure  block  (N)  and  FD  = 
resultant  driving  force  acting  on  failure  block  (N). 

The  resultant  driving  force  acting  on  the  incipient  failure  block  is  the  resultant  of  the 
component  of  the  weight  of  the  failure  block  directed  down  the  failure  plane  and  the 
resultant  of  the  hydrostatic  confining  pressure  directed  up  the  failure  plane: 

FD  =  sinp  -  Fgp  sina  (2) 

where  =  weight  of  failure  block  (Nm'^),  P  =  failure  plane  angle  (degrees),  F^p  = 
magnitude  of  the  resultant  hydrostatic  confining  pressure  (Nm'  )  and  a  =  angle  between 
direction  of  resultant  of  the  confining  pressure  and  a  normal  to  the  failure  plane  (degrees). 

The  resisting  forces  acting  on  the  failure  block  are  a  function  of  the  cohesive  and 
frictional  resistance  to  sliding  along  the  failure  plane.  Accounting  for  the  effects  of  pore 
water  and  confining  pressures,  the  resisting  force  is  given  by: 


FR  =  c  L  +  ((W^f  -  Uy,)  cosP  +  F(.p  cosa)  tanij) 


(3) 


where  L  =  length  of  failure  plane  (m),  =  total  pore  water  pressure  acting  on  failure 

plane  (Nm'^),  c  =  soil  cohesion  (Nm'^)  and  (])  =  soil  friction  angle  (degrees).  By  geometry: 


\  tanp  ■  tani  ^ 


(4) 


where  y  =  soil  unit  weight  (Nm"^),  H  =  overall  bank  height  (m),  H'  =  uneroded  bank 
height  (m),  K  =  tension  crack  depth  (m),  =  relic  tension  crack  depth  (m)  and  i  =  bank 

angle  (degrees).  The  length  of  the  failure  plane  is  given  by: 


H  -  K 
sinP 


(5) 
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Substituting  (2)  through  (5)  into  (1)  gives: 


FS  =- 


C  (H  -K) 
sinP 


2 


tanp 


tani 


-  C/vv  cosp  +  F(.p  cosa 


tan(j) 


1 

2 


(H^  -  K 


2  h'^-K^ 


tanP 


tani 


— j  sinp  -  i^cpsina 


To  apply  equation  (6)  it  is,  therefore,  necessary  to  estimate  the  failure  plane  angle 
and  pore  water  and  confining  pressure  terms  for  given  bank  geometry,  soil  property  and 
ground  and  surface  water  characteristics.  Additionally,  it  is  also  necessaiy  to  determine  the 
location  of  the  most  critical  failure  surface,  in  order  to  relax  the  constraint  that  the  failure 
plane  must  pass  through  the  toe  of  the  bank.  Each  of  these  factors  are  now  considered  in 
turn. 


3.2  Pore  Water  Pressure 

In  general,  the  pore  water  pressure  exerted  on  the  incipient  failure  plane  at  any  point 
is  given  by: 

»<w  =  Pw  g  (K  +  0  (7) 

where  PO^e  water  pressure  at  a  point  on  the  failure  plane  (Nm'^),  =  density  of 

•2 

water  (kgm‘^),  g  =  gravitational  acceleration  (9.81  ms  ),  =  groundwater  head  (m)  and  q 

=  surface  water  surcharge  (m)  (Figure  5).  Inclusion  of  the  surface  water  surcharge  term 
accounts  for  the  influence  of  interacting  pore  and  surface  water  pressures  on  the  overall 
pore  pressure  distribution.  So,  the  pore  water  pressure  distribution  along  the  incipient 
failure  plane  may  be  determined  if  the  phreatic  and  surface  water  surfaces  are  known.  To 
determine  the  total  pore  pressure  acting  on  the  length  of  the  incipient  failure  plane,  it  is 
necessary  to  integrate  (7)  along  the  failure  plane  length,  so  that: 

X 

~  ^  (8) 

0 

where  =  total  pore  water  pressure  acting  on  the  failure  plane  (Pa). 

So,  for  any  phreatic  surface  and  surface  water  distributions,  may  be  estimated 
by  expressing  and  q  as  functions  of  the  phreatic  and  surface  water  elevations,  the 
geometry  of  the  incipient  failure  block,  and  x,  the  total  failure  block  width  (m).  This  can  be 
done  numerically  for  arbitrary  distributions  of  ground  and  surface  water  distributions. 
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However,  it  is  also  possible  to  obtain  numerical-analytical  solutions  for  assumed  special 
cases  of  bank  and  phreatic/surface  water  geometries.  This  is  the  approach  adopted  here.  In 
this  analysis,  it  is  assumed  that  the  phreatic  surface  is  horizontal  and  parallel  to  the  flood 
plain  surface.  For  this  case,  it  is  possible  to  formulate  for  a  range  of  assumed  ground 
water  and  surface  water  elevations.  A  total  of  1 1  combinations  are  possible.  These  are 
illustrated  in  Figure  6  which,  for  clarity,  also  provides  the  various  solutions  for  the  total 
pore  water  pressure  acting  on  the  failure  plane. 


Eigure  5  Impact  of  surface  and  ground  water  on  pore  pressure  at  a  point  on  the  incipient 
failure  plane 

It  is  important  to  stress  that  the  solutions  given  in  Figure  6  represent  pore  water 
pressures  acting  along  the  failure  plane  surface  when  the  pore  and  surface  waters  are 
dynamically  interacting.  This  is  a  departure  from  analyses  which  utilise  pore  pressure  ratios 
to  characterise  the  impact  of  pore  water  pressures  on  bank  stability  (e.g.  Huang,  1983; 
Simon  et  aL,  1991).  Such  analyses  are  not  able  to  account  for  the  impacts  of  surface  water 
surcharge  on  pore  water  pressure.  It  should  be  noted  that  the  effects  on  bank  stability  of 
hydrostatic  thrust  exerted  on  a  tension  crack  filled  by  precipitation  are  not  considered  here. 
Tension  cracks  filled  by  groundwater  have  no  net  hydrostatic  thrusts  exerted  on  them. 


3.3  Hydrostatic  Confining  Pressure 

Submerged  incipient  failure  bocks  may  be  subjected  to  a  hydrostatic  pressure 
exerted  by  the  surface  water  in  the  channel.  The  aim  of  this  section  is  to  determine  the 


SOLUTION  6A  SOLUTION  6D 


SOLUTION  6B 


SOLUTION  6E 


SOLUTION  6G 

(yk<GWSE.^y,p;y3  <y^syj;  GWSE  >y,) 
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SOLUTION  6H 


s 


hydrostatic  confining  pressure  acting  on  submerged  failure  blocks  with  the  geometry 
shown  in  Figure  3  for  arbitrary  surface  water  and  phreatic  surface  elevations.  A  similar 
procedure  was  used  by  Simon  et  al.  (1991),  but  they  considered  banks  only  with  the 
simple  geometry  shown  in  Figure  2.  They  also  did  not  account  for  interactions  between  the 
pore  and  surface  waters  in  estimating  the  resultant  confining  pressure.  These  limitations  are 
addressed  here. 


Figure  7  Influence  of  hydrostatie  confining  pressure  on  bank  stability  in  absence  of 
ground  water 

Consider  Figure  7,  in  which  the  influence  of  the  phreatic  surface  is  neglected.  The 
parameters  of  interest  from  the  point  of  view  of  impact  on  bank  stability  (see  Figure  4)  are 
the  magnitude  of  the  resultant  confining  pressure,  F(.p,  and  the  angle,  co  (degrees),  at  which 
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the  resultant  is  directed  through  the  bank  surface,  since  these  determine  the  components  of 
Fgp  acting  on  the  incipient  failure  block,  via  the  angle  a.  In  fact  the  resultant  hydrostatic 
confining  pressure  acting  on  an  immersed,  inclined,  plane  with  complex  geometry,  such  as 
that  shown  in  Figure  7,  can  be  estimated  by  resolving  the  forces  acting  horizontally  and 
vertically  on  the  immersed  bank  surface: 

fcp  =  VF*^  +  Fy^  (9) 


The  angle  at  which  this  resultant  is  directed  through  the  bank  surface  is  given  by: 


tan  CO  =  ^ 


(10) 


where  =  horizontal  component  of  confining  pressure  (Nm'^)  and  Fy  =  vertical 
component  of  confining  pressure  (Nm‘^).  The  total  horizontal  and  vertical  forces  are  given 
by  the  sum  of  all  the  forces  acting  on  small,  sub-elements  of  the  bank  surface: 


(11) 

Fy  =  13Fy  =  W 

(12) 

where  W  =  weight  of  the  water  above  the  bank  surface  (Nm'^),  that  is  the  weight  of  the 
water  contained  in  the  volume  BCDE  in  Figure  7.  Fy  is,  therefore,  given  by: 

Fy  =  Vk  =  g  (Area  BCDE)  (13) 

The  precise  value  of  Area  BCDE  depends  on  the  surface  water  elevation,  so  the  following 
possibilities  are  obtained: 


Fy  =  ^  cot/  [(y^  -  ys)2  -  (y^  - 

()'w>>'t) 

(14) 

Fy=^cot/(yw-ys)^ 

Cys<3’w^yt) 

(15) 

Fy  =  0 

CO 

VI 

(16) 

where  =  unit  weight  of  the  surface  water  (assumed  to  be  9810  Nm'^)  and  yf,  and  = 
elevations  of  points  on  the  bank  profile  (m)  (see  Figure  7)  and  y^  =  water  surface  elevation 
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(m).  The  value  of  is  given  by  the  sum  of  all  horizontal  components  acting  on  the  bank 
face,  so  that: 

^X  =  '^xl+^x2+^x3  (17) 

r\ 

where  =  pressure  force  acting  on  projection  of  face  AC  (Nm'  ),  Fx2  =  pressure  force 

,  o 

acting  on  projection  of  face  CD  (Nm  )  and  Fx3  =  pressure  force  acting  on  projection  of 
face  DF  (Nm  ).  These  terms  are  given  by  the  following  set  of  equations: 


7^xi=^(yw-3’t)^ 

(yw>>’t) 

(18) 

0 

II 

(>’w^3’t) 

(19) 

P:a=‘\<y^-ys9 

(y^^yt) 

(20) 

<y^>yt) 

(21) 

F:a=\<ys-yd^ 

Al 

(22) 

^x3=^0’w-yt)^ 

(yw<ys) 

(23) 

these  expressions  can  be  substituted  back  into  equations  to  obtain  the  following 
expressions: 


=  ^/(^  (^’w  -  coti  [(y^  -  -  (y^  -  y^)^] 


tan  © 


coti  [(y^  -  ys)  -  (yw  -  >^t)  1 
(^w  - 


CVw  >  3^1) 


(24) 


(25) 


(26) 
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(J-W  -  yi)^ 

VI 

VI 

(27) 

^cp=^  (yw-w)^ 

(28) 

tan  co  =  0 

(>'w<Js) 

(29) 

o 

II 

(30) 

tan  CO  =  0 

VI 

(31) 

The  angle,  O),  at  which  the  resultant  of  the  confining  pressure  passes  through  the  bank 
surface  controls  the  angle,  a  (degrees),  between  the  resultant  of  the  confining  pressure  and 
the  normal  to  the  incipient  failure  plane  (Figure  4),  which  influences  the  components  of 

acting  on  the  incipient  failure  block.  By  geometry,  this  angle  is  given  by: 

a  =  90  -  (P  +  CO) 

(32A) 

a  =  (/  -  P) 

(32B) 

where  equation  (32B)  is  used  only  for  submerged  bank  surfaces  which  are  planar  in  shape 
(i.e.  H  =  H'\  =  0). 

It  is  possible  to  show  that  the  line  of  action  of  the  resultant  of  the  confining  pressure 
never  intersects  the  tension  crack  (Figure  4),  so  long  as  the  following  condition  is  satisfied: 

H 

H  -  K  >-j - X  tan  Cd  (33) 

This  condition  is  not  likely  to  be  violated,  since  tension  cracks  depths  rarely  exceed  half  the 
overall  bank  height  (Taylor,  1948).  This  is  not  a  trivial  result,  for  this  guarantees  that  the 
line  of  action  of  the  resultant  of  the  confining  pressure  has  the  impact  on  the  force 
distribution  on  the  incipient  failure  block  described  by  the  terms  in  the  factor  of  safety 
equation  (6).  No  modification  of  this  equation  is,  therefore,  required  to  account  for  the 
influence  of  confining  pressure  on  tension  cracking. 

The  equations  developed  above  apply  only  to  the  case  when  there  is  no  phreatic 
surface  present,  or  the  phreatic  surface  is  below  the  base  of  the  incipient  failure  block.  In  all 
other  instances,  the  pore  and  surface  waters  along  the  saturated  portion  of  the  bank  surface 
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interact.  Along  this  portion  of  the  bank  surface  the  hydrostatic  pressure  is  exerted  equally  in 
all  directions,  so  that  the  net  resultant  hydrostatic  confining  pressure  is  equal  to  zero.  The 
resultant  of  the  confining  pressure  must,  therefore,  be  calculated  only  along  the  section  of 
the  submerged  bank  face  between  the  levels  of  the  surface  and  groundwater  elevations 
(Figure  8).  However,  it  is  relatively  simple  to  estimate  the  resultant  hydrostatic  confining 
pressures  under  these  circumstances.  All  that  is  required  is  to  replace  yf  in  equations  (24  to 
31)  with  the  elevation  of  the  phreatic  surface  and  to  revise  the  effective  locations  of  yg  and 
accordingly.  In  fact  there  are  a  total  of  12  possible  combinations  of  surface  water 
elevation,  bank  geometry  and  phreatic  surface  elevation,  with  specific  solutions  for  F(.p,  co 
and  a  in  each  case.  These  cases  are,  for  clarity,  summarized  in  Figure  8. 


SOLUTION  8A  SOLUTION  8D 


SOLUTION  8B 


GWSEsy,:y,  >WSEjy5:y,< 


SOLUTION  8C 


SOLUnON  8E 


SOLUTION  8F 


Figure  8  Solutions  for  hydrostatic  confining  pressure 
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SOLUTION  80 


SOLUTION  8J 


SOLUTION  81 


SOLUTION  8L 


Solutions  for  hydrostatic  confining  pressure 


In  obtaining  solutions  for  the  confining  pressure  terms  it  is  necessary  to  recognise 
the  assumptions  used  in  deriving  the  equations  shown  in  Figure  8.  The  phreatic  and  water 
surfaces  are  assumed  to  be  horizontal,  and  sttictly  speaking  the  equations  are  valid  only  for 
surface  water  with  hydrostatic  pressure  distributions.  Although  this  is  a  good 
approximation,  even  in  many  fluids  undergoing  motion,  it  should  be  recognised  that 
pressure  distributions  close  to  the  riverbanks  may  not  always  be  hydrostatic,  particularly  in 
reaches  with  strong  near  bank  secondary  circulations  with  large  vertical  velocity 
components,  as  may  be  found  close  to  the  banks  in  curved  reaches  or  at  high  flood  flows 
(Darby  and  Thome,  1993) 


3.4  Failure  Plane  Angle 

The  failure  plane  angle,  P,  corresponds  to  that  plane  on  which  the  cohesion,  c,  is 
fully  developed  at  critical  conditions  (Taylor,  1948;  Spangler  and  Handy,  1982;  Osman  and 
Thome,  1988).  This  may  be  determined  by  equating  the  first  derivative  of  the  cohesion 
with  respect  to  the  failure  plane  angle  in  equation  (6)  to  zero.  In  other  words,  the  failure 
plane  angle  may  be  determined  by  estimating  the  angle,  P,  that  satisfies  the  condition: 


ap 


=  0 


(34) 


Rearranging  (6)  in  terms  of  c  gives  the  expression: 

(VFt  sin  p  -  Fgp  sina.sinP)-[(Vft  -  l/^,)cosp.sinp.tan(|)  +  cosa.sinp.tan(|)] 

H- K 


where  W^,  F^p,  a,  and  C/^  are  as  defined  previously.  In  fact,  equation  (35)  is  a  more 
complex  function  of  failure  plane  angle  even  than  it  first  appears,  since  a  and  are 
themselves  functions  of  p.  It  is,  therefore,  necessary  to  rearrange  the  expressions  for  these 
parameters  obtained  in  the  previous  sections  in  terms  of  P  prior  to  estimating  the  first 
derivative  of  (35)  with  respect  to  P  in  order  to  obtain  the  condition  defined  by  (34) 

Differentiating  the  rearranged  version  of  (35)  and  equating  the  result  to  zero  leads  to 
the  following  general  expression,  satisfied  only  by  the  tme  value  of  P: 


ap 


(H  -  K)  (-sin^P  +  cos^ 


+ 


(//'^  -  cosp.sinp^A 

(36) 


26 


where: 


X  =  -cosp.sin  (-90  +  P  -i-  (o)  -  sinp.cos  (-90  -i-  P  +  co)  (37 A) 

X  =  -cos  (i  -  P).sinp  +  sin  (/  -  P).cosP  (37B) 

Y  =  cosp.cos  (-90  4-  P  +  CO)  -  sinp.sin  (-90  -i-  P  +  co)  (38A) 

Y  =  cosp.cos  (/  -  p)  +  sinp.sin  (i  -  P)  (38B) 


where  equations  (37B)  and  (38B)  are  used  only  if  equation  (32B)  was  selected,  and: 


0((W^t  ~  ^w)  cosP.sinp.tan(|)) 


H-  K 


It  is  apparent  that  the  form  of  equation  (36)  depends  on  the  pore  water  and 
hydrostatic  confining  pressure  distributions  acting  on  the  incipient  failure  block.  It  is 
possible  to  show  that,  for  the  special  case  when  H  =  H'  and  when  K  =  Kl^  =  F(>p  =  =  0, 

equation  (36)  yields: 

P - 2^  (40) 


This  is  in  fact  the  basis  for  the  common  usage  of  this  simplified  equation  in  many  previous 
analyses  (e.g.  Lohnes  and  Handy,  1968;  Huang,  1983;  Simon  et  al.,  1991).  However,  it 
will  be  clear  that  the  combination  of  conditions  that  lead  to  the  derivation  of  (40)  will  rarely 
be  satisfied  for  natural  riverbanks,  so  that  analyses  which  employ  (40)  will  have  limited 
validity.  Osman  and  Thome  (1988)  solved  for  the  failure  plane  angle  using  an  approach 
based  on  (34),  but  they  did  not  include  pore  water  and  hydrostatic  confining  pressure  terms 
in  their  analysis,  either  for  factor  of  safety  or  failure  plane  angle.  The  analysis  developed 
herein,  therefore,  represents  a  much  more  complete  and  physically-based  approach  to 
prediction  of  failure  plane  angle  than  previous  analyses. 

To  solve  equation  (36)  for  the  failure  plane  angle  requires  an  iterative  solution, 
since  it  is  necessary  to  formulate  F^p  and  in  terms  of  the  failure  plane  angle,  resulting  in 
a  complex  set  of  expressions.  The  approach  adopted  here  is  to  use  the  Newton-Raphson 
iteration  scheme  in  which  successive  estimates  of  (3  are  obtained  using: 


P  =  Pi-1  - 


F(Pi-i) 

F’(Pi.i) 


(41) 
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where  F(Pi.i)  is  given  by  equation  (36)  and  the  derivative  of  this  (with  respect  to  P)  is 
given  by: 


''(Pi-i)  =  ((-2y  {H  -  K)  (sinp.cosP))  -  (- 
X'\  (_  .  Fen  tan({>  y 


'Y(//- 


JUIL 

H  -  K 


H  -  K 


KY^^)  (-sin'^P  +  cos^P) 
tani.//  -  K 

(42) 


where: 


X'  =  2.sinp.cos  (-90  +  P  +  co)  -  2.cosp.cos  (-90  -t-  P  +  co) 

X'  =  -2  sinp.sin  (i  -  P)  -  2  cosP.cos  (i  -  P) 

Y  =  -2.cosp.sin  (-90  +  P  co)  -  2.sinp.cos  (-90  +  p  -i-  co) 

Y  =  2  cosp.sin  (/  -  P)  -  2  sinp.cos  (i  -  P) 

where  equations  (43B)  and  (44B)  are  used  only  if  equation  (32B)  is  selected  and: 

t  0^((l^t  ~  ^w)  cosp.sinp.tan(|))  1 
""  ap2  H-K 


(43A) 

(43B) 

(44A) 

(44B) 


(45) 


An  initial  estimate  of  the  failure  plane  angle  is  made  using  the  approximation: 


B 

Po  2 


(46) 


Convergence  is  obtained  when  equation  (36)  is  satisfied  within  a  small  error  tolerance 
(0.01).  The  Newton-Raphson  algorithm  is  a  robust  numerical  methexi,  and  convergence  is 
obtained  satisfactorily  in  nearly  all  cases.  However,  certain  combinations  of  F^p,  and 
bank  geometry  lead  to  instances  where  converge  does  not  occur,  or  the  scheme  converges 
to  values  of  failure  plane  angle  which  are  either  negative  or  which  are  greater  than  the  bank 
angle.  If  the  iteration  scheme  fails  to  converge  to  a  realistic  solution,  failure  plane  angle  is 
estimated  using  equation  (46).  Fortunately,  in  test  cases  undertaken  as  part  of  this  study 
such  instances  have  been  found  to  be  very  rare  and  have  been  confined  to  very  stable  banks 
with  failure  surfaces  at  high  elevations,  and  very  large  applied  hydrostatic  confining 
pressures.  The  fact  that  these  convergence  problems  arise  only  on  the  most  stable  failure 
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surfaces,  with  factors  of  safety  well  above  unity  (zero  probability  of  failure)  means  that  this 
unavoidable  numerical  problem  does  not  influence  prediction  of  critical  failure  surfaces. 

There  are  11  possible  solutions  for  the  above  equations  since,  as  discussed 
previously,  there  are  1 1  possible  solutions  of  the  pore  water  pressure  term,  each  of  which 
is  a  distinct  function  of  p.  These  solutions  are  not  presented  here  for  the  sake  of  clarity, 
though  each  possibility  is  written  into  the  source  code  of  the  computer  model  based  on  this 
analysis  (Appendix  IV),  In  any  case,  each  of  the  parameters  in  these  equations  have  been 
presented  in  other  forms  in  the  preceding  sections. 

3.5  Location  of  Most  Critical  Failure  Surface 

The  location  of  the  most  critical  failure  surface  is  predicted  using  an  iterative  search 
scheme  based  on  that  presented  by  Simon  et  al.  (1991).  They  stress  that  this  modification 
represents  an  important  departure  from  conventional  analyses  of  planar  failures.  In 
particular,  this  approach  relaxes  the  assumption  used  in  previous  analyses  {e.g.  Lohnes  and 
Handy,  1968;  Osman  and  Thome,  1988)  that  the  failure  plane  must  pass  through  the  toe  of 
the  bank.  This  is  important,  as  there  is  ample  empirical  evidence  to  indicate  that  mass 
failures  can  occur  well  above  the  bank  toe  (Thome  et  al,  1981;  Simon  and  Hupp,  1986; 
Appendix  HI  this  report). 

The  basis  of  the  computational  scheme  involves  repeated  calculation  of  the  factor  of 
safety  at  a  number  of  different  locations  of  the  intersection  of  the  failure  plane  with  the  bank 
surface.  The  most  critical  failure  surface  is  discriminated  by  identifying  the  minimum 
predicted  factor  of  safety.  The  computational  procedure  is  illustrated  in  Figure  9.  First,  the 
search  procedure  divides  the  bank  profile  into  20  equal  increments  of  bank  height  and 
projects  these  to  the  bank  surface.  Then,  for  each  of  these  points  the  dimensions  of  the 
incipient  failure  block  are  determined  by  projecting  a  potential  failure  plane  (at  angle  (3, 
determined  using  the  methods  described  above)  to  the  flood  plain  surface.  It  is  important  to 
stress  that,  unlike  the  analysis  adopted  by  Simon  et  al  (1991),  P  is  allowed  to  vary  and  is 
calculated  separately  for  each  projected  failure  plane.  Pore  water  and  hydrostatic  confining 
pressures  and  the  factor  of  safety  are  then  all  estimated  using  the  procedures  presented  in 
the  preceding  sections.  The  minimum  factor  of  safety  for  the  initial  20  points  is  then 
located. 

If  the  minimum  predicted  factor  of  safety  is  greater  than  unity,  or  maximum 
predicted  probability  of  failure  (see  below)  is  equal  to  zero  it  is  assumed  that  the  bank  is 
stable  with  respect  to  mass  failure.  Otherwise  a  second  search  procedure  is  used  to  locate 
the  most  critical  surface  with  greater  precision  (Figure  9).  In  this  second  search,  the  vertical 
bank  height  between  the  points  immediately  above  and  below  the  point  at  which  the  most 
critical  surface  located  in  the  first  search  was  identified,  is  divided  into  a  further  20  equal 
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sub-increments.  If  the  most  critical  failure  plane  identified  in  the  first  search  is  located  at  the 
bank  toe,  then  the  original  increment  between  the  bank  toe  and  the  first  point  above  the  toe 
is  divided  into  20  sub-increments.  Computations  are  then  repeated  with  the  most  critical 
failure  surface  located  at  the  sub-point  corresponding  to  minimum  calculated  factor  of 
safety. 


START 


STOP 


laranCTa 


3.6  Geometries  of  Mass  Failures 


If  the  factor  of  safety  falls  below  the  critical  value,  failure  is  predicted  to  occur.  In 
calculating  the  geometry  of  the  failure  block,  it  is  necessary  to  revise  the  effective  bank 
geometry  characteristics  according  to  the  location  of  the  most  critical  failure  surface,  so  that; 
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(47) 


H  =  y{p-yf 
^'=3’fp-3's 
^’=yfp-3'f 
^=}’fp-yk 


(yf<>'s) 

(48) 

(yf^^s) 

(49) 

(yf<yk) 

(50) 

Al 

<4-1 

(51) 

In  modelling  channel  evolution  through  coupling  bank  stability  analyses  with 
numerical  aggradation-degradation  models,  estimates  of  the  total  volumetric  inflow  of 
sediment  from  bank  erosion  along  the  length  of  modelled  reaches  are  required  as  a  loading 
term  in  the  sediment  continuity  equation  (Darby,  1994).  However,  2-dimensional  bank 
stability  theories  of  the  type  described  above  provide  estimates  only  of  failure  block 
dimensions  per  unit  reach  length.  No  information  about  the  longitudinal  extent  of  bank 
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failures  is  provided  by  these  analyses.  In  the  absence  of  information  about  the  longitudinal 
extent  of  mass-failures  the  loading  term  in  the  sediment  continuity  equation  is  estimated 
using; 


FAX 

AT 


(57) 


where  =  total  volumetric  bank  material  inflow  flux  due  to  mass  failure  (m^s'^),  Ax  = 
length  of  model  reach  (m)  and  AT  =  chosen  time  step  over  which  the  total  volumetric 
inflow  of  bank  materials  is  time  averaged  (s).  Application  of  equation  (57)  assumes  that 
once  failure  is  detected  at  the  computational  node  representing  the  model  reach,  banks  along 
the  entire  length  of  the  reach  fail  simultaneously.  However,  single  mass  failures  are  rarely 
observed  over  bank  lengths  of  more  than  a  few  metres  or  tens  of  metres,  even  along  the 
world's  largest  rivers.  The  spatial  extent  of  natural  mass  failures  will,  therefore,  be  at  least 
an  order  of  magnitude  smaller  than  the  scale  of  mass  failures  modelled  using  this  approach. 
Consequently,  the  bank  material  inflow  rate  determined  using  equation  (57)  will  represent 
an  unrealistically  large  maximum  possible  value  which  may  be  termed  the  ''potential"  bank 
material  flux.  Similarly,  the  flood  plain  widening  increment  and  total  failure  block  width 
values  determined  using  equations  (54)  or  (55)  will  also  represent  maximum  "potential" 
values,  BW*  and  x*. 


Table  1  River  bank  soil  and  geometry  characteristics  at  varying  scales  (from  Darbv.  19941 


SCALE 

LENGTH  (km) 

SOILS 

GEOMETRY 

Micro 

(Single  Failure) 

<0.05 

Homogeneous 
(Low  Variability) 

Uniform 

Meso 

(River  Reach) 

0.1  -  10 

Heterogeneous 

Uniform 

(Low  Variability) 

Macro 

(River  Basin) 

>  10 

Heterogeneous 

Non-Uniform 

The  problem  appears  to  lie  in  mis-application  of  2-dimensional  bank  stability 
theories,  designed  to  predict  the  stability  of  relatively  small  lengths  of  bank  in  which  the 
geometrical  and  geotechnical  properties  controlling  the  stability  of  the  bank  can  reasonably 
be  considered  uniform  (Table  1).  Determination  of  the  reach-scale  sediment  loading  term 
requires  application  of  the  bank  stability  theory  over  a  much  larger  spatial  extent  in  which 
the  geotechnical,  if  not  the  geometrical,  properties  cannot  be  expected  to  be  uniform  (Table 
1).  It  follows  that  it  is  possible  to  obtain  more  realistic  predictions  of  reach-scale  stability 


with  respect  to  mass  failure  and  reach-scale  volumetric  influx  of  bank  materials  by  putting 
the  deterministic  bank  stability  theories  outlined  in  the  preceding  sections  into  a 
probabilistic  framework,  so  that  the  influence  of  statistical  soil  property  variations  along  the 
reach  on  the  factor  of  safety  can  be  taken  into  account.  By  extension,  the  length  of  reach 
that  is  unstable  with  respect  to  mass  failure  may  be  equated  to  the  estimated  probability  of 
failure.  This  approach  allows  the  variability  in  geotechnical  properties,  assumed  to  control 
the  micro-scale  variation  of  stability  along  the  reach,  to  be  coupled  with  the  meso-scale 
geometry  in  order  to  more  realistically  model  the  influence  of  sub-grid  scale  controls  on 
river  bank  stability  and  the  response  of  the  reach  as  a  whole.  The  hypothesis  that  the 
observed  variation  in  bank  stability  along  a  river  reach  is  due  entirely  to  spatial  variation  in 
the  geotechnical  characteristics  of  the  bank  materials  implies  that  it  is  reasonable  to  assume 
that  the  reach-scale  geometry  is  uniform.  This  is  anyway  required  within  the  conceptual 
framework  of  a  numerical  model  in  which  cross-sections  are  chosen  to  characterize  the 
morphology  of  individual  reaches. 

The  procedure  for  calculating  the  probability  of  failure  is  based  on  that  of  Huang 
(1983).  In  essence,  probability  distributions  representing  the  soils  properties  are  substituted 
for  single  valued  soil  property  variables  in  the  factor  of  safety  equations.  The  first  step  is, 
therefore,  to  determine  the  frequency  distributions  for  each  of  the  soil  properties:  cohesion, 
friction  angle  and  unit  weight.  It  is  possible  to  define  these  distributions  by  measurement 
(e.g.  Simon,  1989).  By  dividing  each  continuous  frequency  distribution  into  discrete 
classes,  it  is  possible  to  define  a  finite  number  of  combinations  of  soil  properties,  with 
representative  values  of  each  of  the  soil  properties  for  each  class.  Each  of  these 
combinations,  together  with  the  reach-scale  geometry,  may  be  directly  applied  in  the  bank 
stability  theories  to  determine  the  factor  of  safety  for  that  individual  combination  of  soil 
properties.The  probability  of  occurrence  of  that  predicted  factor  of  safety,  PCF^)  is  equated 
to  the  probability  of  occurrence  of  the  soil  property  combination  that  gives  rise  to  that 
particular  factor  of  safety,  given  the  constraints  of  the  specified  bank  geometry: 


P(F5)p  =  P(c)i  *  m)j  *  P(Y)^  (58) 

where  P(F5)jy^  =  probability  of  the  predicted  factor  of  safety  for  an  individual  soil  property 
combination,  P(c),-  =  probability  of  the  cohesion  being  in  the  class  represented  by  the  value 
Cl,  and  P(<l))y  and  P(y)^  =  probabilities  of  the  friction  angle  and  soil  unit  weight  being 
represented  by  the  values  and  Yjfc,  respectively.  If  the  factor  of  safety  calculated  for  a 
specific  soil  property  combination  is  less  than  unity,  failure  is  predicted  for  that 
combination.  If  failure  is  predicted,  the  probability  given  by  equation  (58)  is  stored  while 
computations  proceed  through  all  the  other  possible  combinations  of  soil  properties.  The 
probability  of  the  overall  factor  of  safety  being  less  than  unity  is  found  from  the  sum  of  all 
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previously  calculated  probabilities  corresponding  to  individual  factors  of  safety  of  less  than 
unity: 


RP(F5'<1)  =  XP(/'<^<1)  (59) 

1=1 

where  RP(F5<1)  =  reach-scale  probability  of  failure  and  P(F5<1)  =  probability  of  an 
individual  factor  of  safety  being  less  than  unity. 

Since  the  probability  of  failure  is  assumed  equal  to  the  fraction  of  the  reach  that 
fails,  the  "potential"  bank  flux  and  flood  plain  retreat  increments  determined  from  equations 
(54  to  56)  can  be  modified  to  give  the  more  realistic  values: 


1/^  =  RP(FS<1)  V 

(60) 

BW  =  RP(F5<1)  BW* 

(61) 

x=  RP(FS<l)x* 

(62) 

9  1 

where  =  total  mass  failure  bank  material  influx  per  unit  length  of  reach  (m  s  )  and  BW 
and  X  =  reach-averaged  flood  plain  widening  increment  and  reach-averaged  total  failure 
block  widths,  respectively  (m). 

Following  failure,  the  bank  geometry  is  re-calculated  using  a  reach-averaged 
updating  scheme,  so  that  a  uniform  reach-scale  geometry  is  maintained  at  each 
computational  node.  For  planar  failures,  bank  geometry  is  updated  using: 


H  =  H 

(63) 

H'  =  (H’  (1-RP(F5'<1))  +  (H  RP(F5<1)) 

(64) 

K  =  (,K  (1-RP(FS<1)) 

(65) 

Kh=(Ki,a-RP(FS<l)) 

(66) 

i  =  Q  (1-RP(F5<1))  +  (p  RP(F5<1)) 

(67) 

The  computational  procedure  used  to  calculate  the  probability  of  failure  and  update  the  bank 
geometry  is  summarized  in  Figure  10.  It  will  be  noted  that  the  implementation  of  the  search 
for  the  most  critical  failure  plane  is  similar  to  that  employed  in  the  non-probabilistic 

34 


approach,  except  that  the  most  critical  failure  plane  corresponds  to  the  failure  plane  with  the 
highest  predicted  probability  of  failure. 


RETURN  TO  MAIN 
PROGRAM 


figure  10  Computational  procedure  used  in  probabilistic  bank  stability  analysis  (after 


Darby.  19941 


The  probabilistic  approach  does  not  require  any  assumptions  to  be  made  concerning 
the  frequency  distributions  of  the  alluvial  soils  properties.  In  fact,  these  distributions  may 
be  determined  by  repeated  measurements  of  soil  properties  in  the  field  over  the  scale  of  the 
river  reach  (e.g.  Simon,  1989).  In  effect,  the  new  approach  takes  account  of  observed 


variations  in  soil  properties  deterministically,  emphasizing  that  the  method  is  not  a 
stochastic  method  of  predicting  the  stability  of  river  banks.  The  probabilistic  approach  may 
also,  therefore,  be  used  to  assess  the  stability  of  heterogeneous  bank  materials,  as  well  as 
the  stability  and  land  loss  and  sediment  yield  along  lengths  of  riverbank. 

The  criterion  for  deciding  to  use  the  2-dimensional  or  probabilistic  analysis  is 
somewhat  arbitrary.  If  an  at-a-site  estimate  is  required,  or  if  the  soils  are  homogeneous,  it 
is  best  to  use  the  2-dimensional  analysis.  Over  larger  space  scales,  the  probabilistic  analysis 
is  preferred.  In  the  computer  model  based  on  the  analyses  presented  herein,  the  user  is  able 
to  select  either  of  the  two  analyses  (Appendix  IV). 

3.8  Computational  Procedure 

The  analysis  developed  in  the  preceding  section  has  been  developed  in  the  form  of  a 
computer  model  which  is  supplied,  with  supporting  documentation,  with  this  report 
(Appendix  IV).  The  computational  procedures  followed  in  applying  these  analyses  in  this 
model  are  now  outlined  and  are  illustrated  in  Figure  9. 

First,  the  necessary  bank  geometry,  ground  and  surface  water  and  soil  property 
data  is  supplied  to  the  model.  The  user  then  selects  a  range  of  options  which  control  the 
analysis.  The  bank  profile  is  then  divided  into  increments  for  the  first  sweep  for  the  most 
critical  failure  plane.  In  re-positioning  the  base  of  the  failure  plane  at  each  increment,  the 
effective  bank  geometry  for  the  failure  plane  is  calculated  using  equations  (47)  to  (53).  The 
confining  pressure  terms  are  then  calculated.  The  model  then  proceeds  to  the  bank  stability 
analyses,  depending  on  which  analysis  the  user  specifies.  In  both  cases,  the  computational 
scheme  is  broadly  the  same.  The  confining  pressure  terms,  together  with  the  bank 
geomeuy  and  soil  property  parameters,  allow  the  failure  plane  angle  to  be  estimated.  The 
pore  water  pressure  term  is  then  estimated.  Next  the  factor  of  safety,  or  probability  of 
failure  (depending  on  the  selected  analysis)  are  predicted.  This  procedure  is  repeated  for  all 
the  selected  failure  plane  locations,  enabling  the  location  of  the  most  critical  failure  plane  to 
be  selected.  If  required,  the  second  search  sweep  identifies  the  precise  location  of  the  most 
critical  failure  plane.  The  geometries  of  the  simulated  mass-failures  (if  detected)  are  then 
calculated.  The  output  data  is  directed  to  files  and  computations  stop. 

The  bank  stability  computer  model  supplied  with  this  report  may  also  be  used  to 

estimate  the  amount  of  fluvial  erosion  required  to  bring  initially  stable  banks  to  the  point  of 

failure.  By  selecting  the  "critical  geometry"  analysis  option  of  the  computer  program,  the 

program  automatically  supplies  increments  of  lateral  fluvial  erosion  and  bed  degradation  to 

the  initially  specified  bank,  until  the  point  of  failure  is  reached.  The  increments  of  lateral 

erosion  and  bed  degradation  applied  to  the  bank  are  supplied  in  steps  according  to  a  user- 
dW 

specified  ratio,  where  dW  is  the  increment  of  lateral  fluvial  (particle-by-particle) 
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erosion  (m)  and  dW  is  the  increment  of  bed  degradation  at  the  toe  of  the  bank  (m).  In  each 

dW 

step,  0.05m  of  degradation  is  automatically  applied  to  the  bank  profile.  The  parameter 

therefore,  determines  the  magnitude  of  the  of  lateral  erosion  in  each  step.  During  the  bed 
degradation,  user-supplied  water  surface  and  ground  water  surface  elevations  are  also 
lowered  at  appropriate  rates. 

The  results  file  provides  the  user  with  the  amount  of  degradation  required  to 
destabilize  the  banks,  as  well  as  estimates  of  the  critical  bank  height  and  critical  failure 
block  magnitudes.  This  analysis  will,  therefore,  be  particularly  useful  for  analysis  of  bank 
erosion  and  lateral  channel  stability  in  degrading  systems.  For  example,  the  analysis  could 
be  used  as  a  design  tool  to  estimate  elevations  of  grade  control  structures  required  to 
prevent  knickpoint  migration  destabilising  riverbanks  and  generating  excessive  channel 
widening. 

3.9  Limitations  of  Analysis 

It  is  appropriate  to  state  the  limitations  of  the  new  analysis.  The  analysis  is 
applicable  to  steep,  cohesive  riverbanks  which  fail  along  planar  surfaces.  Banks  are 
assumed  to  be  non-layered  and  relatively  homogeneous,  though  heterogeneous  bank 
materials  may  be  accounted  for  using  the  probabilistic  bank  stability  analysis.  The  effects  of 
vegetation  and  seepage  effects  are  excluded.  In  calculating  the  soil  pore  water  pressure,  the 
phreatic  surface  is  assumed  to  be  horizontal  and  parallel  to  the  flood  plain  surface.  Surface 
water  pressure  distributions  adjacent  to  the  bank  profile  are  assumed  to  be  hydrostatic.  The 
influence  of  hydrostatic  thrast  of  (precipitation)  water  filled  tension  cracks  is  also  excluded. 
Tension  cracks  are  limited  to  maximum  depths  of  no  more  than  half  the  height  of  the  bank. 

It  is  important  to  take  heed  of  these  limitations  during  application  of  the  computer 
program  based  on  the  riverbank  stability  analysis  (Appendix  IV),  to  ensure  that 
computational  instabilities  do  not  arise.  Known  problems  are  now  listed: 

1.  When  using  the  probabilistic  analysis,  assigning  a  soil  unit  weight  frequency  distribution 
with  low  magnitudes  of  unit  weight  with  non-zero  probabilities  of  occurrence  may  lead  to 
computed  failure  block  weights  less  than  the  pore  pressure  exerted  on  the  failure  plane 
when  the  magnitude  of  the  computed  pore  water  pressure  is  also  large.  This  can  in  turn  can 
lead  to  negative  values  of  factor  of  safety,  which  are  difficult  to  analyse  physically. 

2.  Low  values  of  bank  angle  may,  with  large  values  of  H'  and  small  values  of  (3,  lead  to 
computed  failure  block  weights  of  less  than  zero,  which  is  impossible.  The  computer 
model  has  been  successfully  applied  on  banks  with  a  wide  range  of  geometries  and  bank 
angles  as  low  as  50  degrees.  This  value  may  be  used  as  a  guideline  for  the  lower  limit  of 
bank  angle  that  may  be  used  in  this  analysis. 
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3.  Values  of  tension  crack  depth  greater  than  half  the  total  bank  height  may  lead  to  negative 
values  of  failure  block  weight,  which  is  impossible.  The  computer  program  based  on  this 
analysis  automatically  constrains  tension  cracks  so  they  do  not  exceed  this  maximum  depth. 
This  should  be  considered  when  interpreting  analyses  for  which  deeper  cracks  were 
specified  by  the  user. 

4.  ANALYSIS  OF  PREDICTIVE  ABILITY  OF  SELECTED 
RIVERBANK  STABILITY  ANALYSES 

In  this  section,  data  obtained  from  natural  stream  channels  undergoing  erosion 
through  mass-failures  are  used  to  test  the  predictive  abilities  of  some  selected  riverbank 
stability  analyses,  including  validation  of  the  analysis  developed  in  the  preceding  section. 
The  comparison  of  the  predictive  abilities  of  the  various  analyses  is  used  to  elucidate  which 
of  the  analyses  is  most  suited  to  simulating  the  stability  and  geometry  of  cohesive  stream 
banks  undergoing  planar  mass-failures.  The  analysis  should,  therefore,  be  helpful  in  two 
regards.  First,  the  analysis  will  aid  river  engineers  in  their  selection  of  bank  stability 
analyses  for  assessing  channel  stability.  Second,  the  analysis  should  assist  in  selecting  the 
most  suitable  mass-wasting  algorithm  for  inclusion  in  numerical  mobile-bed  models,  such 
as  HEC-6,  for  development  of  mobile-width  numerical  modeling  techniques. 

4.1  Selection  of  Bank  Stability  Analyses 

Even  considering  only  the  planar  mechanism  of  mass-failure,  a  large  number  of 
bank  stability  analyses  are  available  for  consideration  when  selecting  analyses  to  subject  to 
assessment  of  predictive  ability.  At  the  outset  of  this  study,  it  was  planned  to  assess  the 
predictive  abilities  of  a  wide  range  of  analyses,  including  those  developed  by  commercial 
software  companies,  academics,  and  other  professionals.  In  this  way,  it  was  hoped  that 
analyses  in  use  by  a  whole  spectrum  of  professionals  would  be  subjected  to  an  independent 
scrutiny  of  their  predictive  abilities.  To  the  knowledge  of  the  authors,  no  such  independent 
scrutiny  has  been  attempted  previously. 

However,  this  proved  to  be  an  overly  optimistic  goal.  Despite  contacting  many 
individuals,  academics,  commercial  software  companies  and  private  consultants,  response 
to  the  authors  requests  to  provide  these  analyses  was  generally  poor.  Many  were  not 
prepared  to  freely  provide  their  codes,  even  solely  for  the  purposes  of  this  research.  The 
resources  of  this  project  were  not  large  enough  to  purchase  these  codes.  The  selected 
analyses  are,  therefore,  restricted  to  stability  analyses  which  are  in  the  public  domain  and 
have  been  fully  documented  in  the  literature.  Nevertheless,  the  results  of  the  predictive 
comparison  of  the  selected  analyses  will  still  be  very  useful,  because  the  analyses  tested  are 
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by  default  well  known  and  widely  used.  Furthermore,  many  commercially  available  codes 
which  were  not  available  for  testing  in  this  analysis  are  anyway  based  on  analyses  similar 
to  those  of  Lohnes  and  Handy  (1968)  and  Huang  (1983),  which  are  included  in  this 
predictive  comparison. 

In  selecting  bank  stability  analyses  for  assessment,  a  number  of  criteria  were  set. 
Selected  analyses  should  be  well  known  and  widely  used.  It  was  hoped  to  obtain  a  range  of 
complexity  in  the  selected  analyses,  in  order  to  allow  assessment  of  whether  or  not 
increasing  the  physical  basis  and  complexity  of  bank  stability  analyses  is  justified  in  terms 
of  improved  predictive  ability.  All  the  selected  analyses  should  be  applicable  to  planar 
failures  in  cohesive  riverbanks.  The  analyses  should  be  directly  applicable  to  riverbank 
environments. 

A  total  of  4  analyses  were  selected  on  this  basis.  These  are  the  analyses  of  Lohnes 
and  Handy  (1968)  (abbreviated  to  LH68),  Huang  (1983)  (H83),  Osman  and  Thome 
(1988)  (OT88)  and  the  analysis  developed  in  this  study  (Darby-Thome)  (DT94).  The  main 
features  of  these  analyses  have  previously  been  described  in  the  introduction,  but  are  also 
summarized  in  Table  2.  It  can  be  seen  that  the  listed  order  of  these  analyses  roughly 
corresponds  to  an  increasing  level  of  complexity.  It  was  hypothesized  that  the  predictive 
abilities  of  the  analyses  would  improve  down  the  listed  order. 

Table  2  Summary  of  characteristics  of  bank  stability  analyses  selected  for  predictive 
comparison 


Analysis 

Geometry 

Failure 

Plane 

Analysis 

Confining 

Pressure 

Pore  Water 

Pressure 

Tension 

Crack 

LH68 

Simple 

■©- 

+ 

II 

CQ. 

No 

No 

Yes 

H83 

Simple 

No 

Yes 

Yes 

OT88 

"Natural" 

Physically- 

Based 

No 

No 

Yes 

DT94 

"Natural" 

Physically- 

Based 

Yes 

Yes 

Yes 

4.2  Collation  of  Data 


Few  efforts  to  test  rigorously  the  predictive  ability  of  bank  stability  analyses  have 
previously  been  made.  In  no  small  measure  this  is  due  to  the  difficulty  in  obtaining  reliable 


measurements  of  the  critical  geometry  of  mass-failures  at  the  point  of  failure.  Allied  to  the 
difficulty  in  obtaining  the  geometries  of  incipient  failure  blocks  lies  the  problem  of 
measuring  the  geotechnical  soil  properties  of  critical  banks. 

However,  in  the  course  of  this  research  a  reliable  data  set  of  51  critical  stream  banks 
with  corresponding  geotechnical  properties  has  now  been  compiled.  This  set  of  data 
represents  a  resource  of  considerable  value.  This  data  set  has  been  provided  in  spreadsheet 
form  on  diskette,  and  is  supported  by  Appendbc  in  of  this  report 

Geometrical  data  relating  to  critical  failure  blocks  were  obtained  in  the  following 
way.  Repeated  cross-section  surveys  along  destabilized  streams  widening  through  bank 
failures  have  for  a  number  of  years  now  been  made  by  United  States  Department  of 
Agriculture  personnel  at  frequent  intervals  along  three  streams  in  northern  Mississippi: 
Hotophia  Creek,  Long  Creek  and  Goodwin  Creek.  By  overlaying  plots  of  the  sequences  of 
the  cross-section  profiles  obtained  from  the  surveys,  it  was  possible  to  estimate  changes  in 
bank  geometry  between  surveys.  In  many  instances  these  changes  in  cross-section 
geometry  were  characterised  by  mass  failures  (Figure  1 1).  In  this  way,  the  geometries  of 
mass  failures  were  reconstructed  using  the  survey  data.  When  inspecting  the  over-laid 
profiles,  care  was  taken  to  ensure  that  changes  in  geometry  observed  between  surveys 
reflected  only  changes  caused  by  mass  failures.  To  ensure  quality  control,  data  were 
discarded  if  there  was  any  doubt  that  an  observed  change  in  the  bank  profile  did  not 
represent  a  planar  mass  failure. 


Figure  11  Example  of  over-laid  cross-section  surveys  used  to  estimate  critical  geometries 
of  mass-failures 


The  bank  profile  data  so  obtained  were  supplemented  by  measurements  of  the 
geotechnical  properties  of  the  bank  materials.  Direct  measurements  of  soil  properties  have 
been  made  at  sites  along  Lx)ng  and  Goodwin  Creek  by  Thome  et  al.,  1981.  Full  details  of 
the  measurement  techniques  and  data  so  obtained  are  available  in  that  report.  It  is  assumed 
here  that  these  measured  properties,  made  only  at  a  few  sites,  reflect  soil  characteristics 
throughout  the  period  in  which  channel  surveys  were  made  (up  to  1988).  Since  the 
stratigraphy  and  soil  units  at  the  Hotophia  Creek  sites  were  generally  similar  to  those  in  the 
adjacent  Long  and  Goodwin  Creek  watersheds,  soil  properties  for  the  Hotophia  Creek  sites 
were  estimated  from  the  values  determined  at  these  nearby  sites  (Thome  et  al,  1981). 

4.3  Assessment  of  Predictive  Ability  of  Stability  Analyses 

Assessment  of  the  predictive  abilities  of  each  of  the  tested  analyses  is  made  in  two 
categories.  First,  the  ability  of  each  method  to  predict  the  overall  stability  of  the  bank  with 
respect  to  mass  failure  under  gravity  is  assessed  by  comparing  predicted  factors  of  safety  to 
the  critical  factor  of  safety  value  of  unity.  Second,  the  ability  of  each  analysis  to  predict  the 
geometry  of  incipient  failure  blocks  is  assessed  by  comparing  predicted  failure  block 
geometries  with  observed  failure  block  geometries.  These  two  aspects  are  now  discussed. 

4.3.1  Factor  of  Safety  Analysis 

Comparisons  of  predicted  and  observed  factor  of  safety  for  each  selected  stability 
analysis  were  made  in  the  following  way.  Observed  critical  bank  geometries  and 
geotechnical  parameters  at  each  of  the  51  sites  in  the  data  set  were  used  as  input  data  in  each 
of  the  4  stability  analyses,  to  obtain  predicted  factors  of  safety  for  each  analysis  at  each  of 
the  51  sites.  These  predicted  values  were  then  compared  to  "observed"  values  of  factor  of 
safety.  For  perfect  agreement,  predicted  factors  of  safety  should,  at  each  site,  be  less  than 
or  equal  to  the  "observed"  critical  value  of  unity.  Assessment  of  the  overall  predictive 
ability  of  each  analysis  can,  therefore,  be  made  by  comparing  discrepancies  between 
predicted  and  observed  factor  of  safety. 

In  comparing  predicted  and  observed  factors  of  safety  in  this  way,  the  reliability  of 
the  assessment  of  the  relative  predictive  abilities  of  each  analysis  is  dependent  upon  the 
quality  of  the  data.  While  it  is  felt  that  the  geometry  of  observed  failures  is  described  with 
high  accuracy  and  reliability  by  this  data  set,  no  direct  observations  of  the  soil  geotechnical 
parameters  were  made  at  the  point  of  failure  at  any  of  the  specific,  individual,  sites  in  the 
data  set.  Instead,  geotechnical  data  were  obtained  from  (reliable)  measurements  of  soil 
properties  at  a  few  sites  representative  of  the  catchments  in  which  the  observed  bank 


failures  were  located.  Moreover,  these  measurements  where  made  when  the  soils  were 
unsaturated  and,  therefore,  represent  "average",  rather  than  the  "worst  case"  conditions 
during  which  mass  failures  most  likely  occurred  (Thome  et  al.,  1981).  "Worst  case" 
estimates  of  the  soil  properties  at  each  site  were  made  by  modifying  arbitrarily  the  measured 
properties.  Predictive  comparisons  were  then  made  for  both  "average"  and  "worst  case" 
conditions.  In  this  way,  although  the  precise  values  of  the  geotechnical  data  at  the  point  of 
failure  are  not  known  for  any  of  the  sites,  the  range  of  values  between  the  estimated  worst 
case  and  average  conditions  almost  certainly  encompass  the  actual  values  encountered. 

Figure  12  illustrates  the  results  of  the  predictive  comparison  for  each  of  the  4 
selected  analyses  for  the  51  sites  combined,  for  both  average  and  worst  case  soil 
properties.  For  comparative  purposes,  the  "observed"  critical  factor  of  safety  (FS  =  1)  is 
also  illustrated.  For  an  analysis  to  have  been  successful  at  any  one  site,  in  terms  of 
predicting  that  the  observed  bank  was  unstable,  the  predicted  factor  of  safety  must  be  less 
than  or  equal  to  this  critical  value.  For  the  purposes  of  this  study,  a  reasonable  tolerance 
level  was  specified,  such  that  predictions  of  factor  of  safety  were  considered  accurate  if  the 
predicted  value  of  factor  of  safety  was  less  than  or  equal  to  1.2.  This  value  is  close  enough 
to  the  critical  case  to  be  considered  as  successfully  predicting  a  bank  of  unreliable  stability. 

(A)  (B) 


Figure  12  Comparison  of  frequencies  of  occurrence  of  predicted  factor  of  safety  and 
"observed"  factor  of  safety  for  (Al  Lohnes  and  Handy  (19681.  (B1  Huang  (19831.  fC 
Osman  and  Thome  (1988).  and  (D)  Darbv-Thome  bank  stability  analyses 


Figure  12  shows  that  all  4  tested  analyses  tend  to  over-predict  factor  of  safety.  That 
is,  all  the  tested  analyses  tend  to  predict  more  stable  banks  than  actually  observed.  All  of 
the  analyses  indicate  correctly  that  factor  of  safety  is  predicted  to  decrease  as  soil  properties 
tend  to  the  "worst  case"  condition.  However,  there  is  considerable  variation  in  both  the 
degree  to  which  each  of  the  selected  analyses  over-predicts  the  factor  of  safety,  and  in  the 
range  of  values  predicted  by  each  analysis  (the  degree  of  scatter  of  predicted  values). 
Assuming  that  the  "worst-case"  data  corresponds  to  the  conditions  when  the  banks  were 
most  likely  to  have  failed  (Thome  et  al.,  1981),  Figure  12  indicates  that  the  most  successful 
analysis  in  terms  of  having  the  highest  number  of  sites  predicted  within  the  required  range 
of  ^5  <  1.2  ("worst-case"  conditions)  is  the  Darby-Thome  analysis  (20  out  of  51  sites  fall 
in  this  category).  The  Darby-Thome  analysis  is  followed  by  the  Osman-Thome  (14  out  of 
the  51  sites),  Lohnes-Handy  (5  sites)  and  Huang  (3  sites)  analyses.  This  ranking  of 
predictive  ability  is  as  expected  and  corresponds  to  the  increasing  level  of  complexity  of  the 
various  analyses  (Table  2). 

To  quantify  an  overall  level  of  predictive  ability  of  each  of  the  bank  stability 
analyses,  comparisons  of  predicted  versus  observed  factors  of  safety  can  be  made  in  terms 
of  the  overall  mean  of  the  ratios  of  predicted  to  "observed"  factors  of  safety  for  each  of  the 
51  sites  in  the  data  set.  This  average  value  provides  a  simple,  but  reliable,  quantitative 
estimate  of  the  overall  difference  between  predicted  and  observed  values  for  each  analysis. 
Assuming  that  the  "observed"  factor  of  safety  is  equal  to  one,  the  average  discrepancies 
obtained  for  each  analysis  represent  the  mean  of  the  51  predicted  factors  of  safety.  These 
values  are  summarised  in  Table  3.  This  table  confirms  the  impression  of  Figure  12  that  the 
Darby-Thome  analysis  is  the  most  successful  of  the  tested  analyses,  but  even  this  analysis 
still  somewhat  over-predicts  factor  of  safety. 


Table  3  Summary  of  predictive  abilities  of  selected  bank  stability  analyses 


Analysis 

Mean  Predicted  Factor 

of  Safety  ("Average" 
soils) 

Mean  Predicted  Factor 

of  Safety  ("Worst- 
case"  soils) 

Darby-Thome  (This  report) 

2.796 

1.431 

Osman  and  Thome  (1988) 

3.200 

1.816 

Lohnes  and  Handy  (1968) 

2.690 

1.828 

Huang  (1983) 

3.676 

3.257 

4.3.2  Geometry  of  Incipient  Failure  Blocks 

Although  the  stability  of  the  banks  with  respect  to  mass  failure  is  important  in 
determining  when  the  bank  will  collapse,  predictions  of  the  failure  block  geometry  must 
also  be  correct  if  the  amount  of  widening,  volume  of  sediments  delivered  to  the  channel  and 
morphology  of  the  bank  following  mass-failure  are  also  to  be  estimated  accurately. 
Knowledge  of  these  parameters  is  essential  if  quantitative  predictions  of  the  magnitudes  of 
channel  adjustments  are  required.  In  this  regard  the  critical  parameters  are  the  failure  plane 
angle,  which  determines  the  angle  of  the  "new"  bank  surface  following  a  mass-failure;  the 
failure  block  volume,  which  determines  the  amount  of  bank  material  delivered  to  the 
channel;  and  the  width  of  the  failure  block,  which  determines  the  magnitude  of  channel 
widening.  The  latter  parameter  is  described  here  using  both  the  total  width  of  the  failure 
block  and  the  flood-plain  widening  increment  (Figure  4). 

Each  of  the  4  bank  stability  analyses  tested  in  this  study  provide  predictions  of  each 
of  the  parameters  described  above.  Indeed  it  is  the  volume  and  shape  of  the  incipient  failure 
block  that  determines  the  diving  and  resisting  forces  which  control  the  stability  of  the  block 
with  respect  to  mass-failure.  Assessment  of  the  ability  of  each  bank  stability  analysis  to 
predict  the  critical  failure  block  geometry,  therefore,  provides  important  clues  regarding  the 
reasons  for  the  success  or  failure  of  the  various  analyses  to  predict  factor  of  safety.  It  is  this 
aspect  that  is  examined  herein.  To  achieve  this  goal,  predicted  incipient  failure  block 
geometries  are  compared  directly  to  the  observed  critical  failure  block  geometries  for  each 
of  the  4  selected  bank  stability  analyses.  Predictive  ability  is  measured  quantitatively  using 
the  mean  of  the  discrepancy  ratio  {Me)  and  the  absolute  deviation  of  the  discrepancy  ratio 
{Ad): 


A/e  =  i5;fE  (68) 

Ad=\-L\^^-Me\\  (69) 

where  Xp  =  value  of  predicted  parameter,  Xq  =  value  of  observed  parameter  and  n  = 
number  of  data  points  in  the  data  set.  It  is  seen  that  the  mean  of  the  discrepancy  ratio  is 
defined  as  the  arithmetic  mean  of  the  ratios  of  predicted  to  observed  parameters  for  the  data 
set  as  a  whole,  while  the  absolute  deviation  of  the  discrepancy  ratio  is  a  measure  of  scatter. 
It  follows  that  perfect  agreement  between  predicted  and  observed  parameters  occurs  when 
Me  =  1  and  Ad  =  0. 

It  should  be  noted  that  this  approach  is  not  a  comparison  of  the  ability  of  the  various 
analyses  to  predict  failure  block  geometry  at  the  point  of  failure.  This  is  because  for  any 
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predicted  failure  block  geometry  obtained  using  the  observed  data  the  predicted  factor  of 
safety  may  not  be  below  the  critical  value,  so  that  failure  is  not  predicted  to  actually  occur 
for  that  incipient  failure  block.  To  compare  the  ability  of  the  models  to  predict  the  critical 
failure  geometry  parameters  it  would  be  necessary  to  use  the  models  to  determine  the 
critical  geometry  and  then  compare  these  predicted  values  with  the  data.  However,  a 
comparison  of  each  of  the  ability  of  each  analysis  to  predict  the  critical  geometry  could  not 
be  attempted  here  because  the  Lohnes-Handy  and  Huang  analyses  are  not  implemented  in 
this  way.  Instead  the  comparison  of  predictive  abilities  of  each  analysis  is  restricted  to  the 
case  described  in  order  to  elucidate  reasons  for  failure  to  predict  accurately  factor  of  safety. 
A  full  analysis  of  the  ability  of  the  Darby-Thome  analysis  to  predict  failure  block  geometry 
at  critical  conditions  is,  however,  provided  in  section  4.3.3. 

Table  4  Summary  of  abilities  of  selected  bank  stability  analyses  to  predict  failure  block 
geometry  ("Worst-Case"  conditions') 


Model 

Failure  Plane 
Angle,  p 

Failure  Block 

Volume,  VB 

Failure  Block 

Width,  X 

— 

Me 

Ad 

Me 

Ad 

Me 

Ad 

Me 

Ad 

LH68 

1.237 

0.386 

0.391 

0.259 

0.703 

0.311 

0 

0 

Huang 

1.281 

1.309 

1.030 

0.736 

1.249 

0.482 

2.053 

1.464 

OT88 

0.561 

0.129 

2.586 

1.417 

1.714 

0.551 

3.478 

1.919 

DT94 

0.867 

0.359 

1.982 

1.372 

1.349 

0.627 

2.951 

1.848 

Table  4  summarises  calculated  discrepancy  ratios  and  absolute  deviation  parameters 
for  the  incipient  failure  block  geometry  parameters  for  each  of  the  tested  bank  stability 
analyses  for  "worst-case"  conditions.  In  addition,  the  scatter  plots  Figures  13  to  15 
illustrate  graphically  the  comparisons  of  predicted  and  observed  values  of  failure  plane 
angle,  failure  block  volume  and  failure  block  width,  again  for  "worst-case"  conditions 
only. 


Predicted  Beta  (Degrees) 


(A)  Lohnes-Handy  Analysis 


(B)  Huang  Analysis 


Observed  Beta  (Degrees) 


Observed  Beta  (Degrees) 


(C)  Osman-Thorne  Analysis 


(D)  Darby-Thorne  Analysis 


Observed  Beta  (Degrees) 


Observed  Beta  (Degrees) 


Figure  13  Comparison  of  predicted  versus  observed  failure  plane  angle  for  (Al  Lohnes- 
Handy.  (B)  Huang.  (C)  Osman-Thorne,  and  (T)')  Darbv-Thome  bank  stability  analyses 
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Predicted  VB  (in^3/in)  Predicted  VB  (in''3/ni) 


(C)  Osman-Thorne  Analysis  (D)  Darby-Thorne  Analysis 


Observed  VB  (in''3/m)  Observed  VB  (m''3/m) 


Figure  14  Comnarison  of  predicted  versus  observed  failure  block  volume  for  ('A'l 
Lohnes-Handv,  (B)  Huang.  (C)  Osman-Thome.  and  (D^  Darhv-Thome  hank  stability 
analyses 
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Figure  13  shows  that  distinct  differences  arise  in  predictions  of  failure  plane  angle 
according  to  the  methods  of  prediction.  The  Osman-Thome  and  Darby-Thome  methods 
both  employ  a  physically-based  procedure  to  predict  failure  plane  angle  with  relatively 
good,  consistent,  results  though  with  a  tendency  to  somewhat  under-predict  failure  plane 
angle.  This  impression  is  confirmed  by  the  calculated  mean  discrepancy  ratios  of  0.561  and 
0.867  for  the  Osman-Thorne  and  Darby-Thorne  analyses,  respectively  (Table  4). 
However,  when  the  simple  approximation  based  on  bank  and  soil  friction  angles  only  is 
used  to  calculate  failure  plane  angle,  as  is  employed  in  both  the  Lohnes-Handy  and  Huang 
analyses,  predictive  ability  decreases  dramatically.  Figure  13  clearly  shows  the  poor 
predictive  ability  of  these  two  analyses  for  failure  plane  angle,  to  an  extent  that  is  masked 
by  the  calculated  mean  discrepancy  ratios  for  these  analyses  (Table  4).  The  mean 
discrepancy  values  of  1.237  (Lohnes-Handy)  1.281  (Huang)  appear  to  suggest  improved 
predictive  ability  relative  to  the  Osman-Thorne  analysis,  but  this  is  misleading.  In  fact 
predictions  for  these  two  analyses  tend  to  be  distributed  equally  badly  either  side  of  the 
correct  values  (Figure  13),  forcing  the  mean  discrepancy  ratio  closer  to  unity  than  is 
actually  warranted.  It  is  important  to  note  that  overall  the  Osman-Thome  and  Darby-Thome 
analyses  tend  to  under-predict  failure  plane  angle,  which  will  tend  to  contribute  to  the  over¬ 
prediction  of  factor  of  safety  observed  in  the  previous  section.  However,  the  Lohnes- 
Handy  and  Huang  analyses  tend  to  over-predict  failure  plane  angle,  suggesting  that  the 
cause  of  over-prediction  of  factor  of  safety  for  these  analyses  must  be  sought  elsewhere.  Of 
the  tested  analyses  the  Darby-Thome  analysis  predicts  incipient  failure  plane  angles  closest 
to  the  observed  data. 

Figure  14  illustrates  that  all  of  the  analyses  tend  to  give  a  large  amount  of  scatter  in 
predicted  versus  observed  failure  block  volume  values.  The  diagrams  also  suggest  that 
while  the  Osman-Thome  and  Darby-Thome  analyses  tend  to  over-predict  failure  block 
volume,  the  Lohnes-Handy  analysis  tends  to  under-predict  this  parameter.  The  Huang 
analysis  appears  to  give  values  which  bisect  the  line  of  perfect  agreement,  but  with  a  high 
degree  of  scatter.  These  qualitative  impressions  are  confirmed  by  the  calculated  values  of 
mean  discrepancy  ratio:  Lohnes-Handy,  0.391;  Huang,  1.030;  Osman-Thome,  2.586;  and 
Darby-Thome,  1.982,  though  the  poor  absolute  deviation  values  (Table  4)  indicate  the  high 
degree  of  scatter  for  all  of  these  analyses.  The  magnitude  of  the  volume  of  the  failure  block 
determines  the  weight  of  the  block  and,  therefore,  the  predicted  factor  of  safety.  Over¬ 
predictions  of  failure  block  volume  for  the  Osman-Thome  and  Darby-Thorne  analyses 
should  tend  to  result  in  under-prediction  of  factor  of  safety  (/.e.  banks  less  stable  than 
observed).  That  the  opposite  is  in  fact  the  case  for  both  these  analyses  indicates  that  under¬ 
prediction  of  failure  plane  angle  compensates  for  over-prediction  of  block  volume  in  terms 
of  its  impact  on  factor  of  safety.  Conversely,  the  over-prediction  of  failure  plane  angle  for 


the  Lohnes-Handy  and  Huang  analyses  is  compensated  by  under-prediction  of  failure  block 
volume  to  the  extent  that  factor  of  safety  is  also  over-predicted  by  both  these  analyses. 
Indeed,  the  failure  plane  angle  also  directly  influences  the  predicted  block  volume;  lower 
angles  leading  to  higher  block  volumes  and  vice  versa.  These  results  suggest  that  prediction 
of  failure  plane  angle  is  critical  in  estimating  the  dimensions  of  the  incipient  failure  block 
and  the  stability  of  the  bank. 

Figure  15  illustrates  the  predicted  versus  observed  failure  block  widths  obtained  for 
each  of  the  four  analyses.  Again,  there  appears  to  be  a  large  amount  of  scatter  in  these 
relationships  for  all  of  the  tested  analyses,  though  there  is  a  hint  that  the  Osman-Thome  and 
Darby-Thome  analyses  tend  to  cluster  in  a  slightly  more  concentrated  form.  The  scatter 
plots  suggest  that  the  Lohnes-Handy  analysis  tends  to  under-predict  failure  block  width, 
while  the  Huang,  Osman-Thome  and  Darby-Thome  analyses  all  appear  to  over-predict  this 
parameter.  This  impression  is  reflected  in  the  calculated  mean  discrepancy  ratios.  In  terms 
of  ranking  predictive  ability  for  this  parameter,  the  Huang  analysis  (1.249)  appears  to  have 
the  mean  discrepancy  ratio  closest  to  unity,  followed,  in  order,  by  Lohnes-Handy  (0.703), 
Darby-Thome  (1.349)  and  Osman-Thome  (1.714).  However,  Figure  15  suggests  that 
there  is  a  tendency  for  the  Darby-Thome  and  Osman-Thome  to  systematically  over-predict 
failure  block  width  which  is  not  apparent  in  either  of  the  other  two  analyses.  The  cause  of 
the  over-prediction  for  these  two  analyses  is  most  likely  the  observed  under-prediction  of 
failure  plane  angle  described  previously. 

Overall,  it  appears  that  the  more  consistent  failure  plane  angle  predictions  of  the 
Darby-Thome  and  Osman-Thome  analyses  enables  these  two  analyses  to  provide  superior 
predictions  of  incipient  failure  block  geometry  and,  therefore,  factor  of  safety  when 
compared  to  the  Lohnes-Handy  and  Huang  analyses.  It  appears  that  the  key  to  improving 
the  predictive  ability  of  the  various  bank  stability  analyses  lies  in  more  accurately  estimating 
failure  plane  angle.  In  this  regard,  it  is  possible  that  the  Darby-Thome  analysis  would 
provide  further  improved  predictions  of  failure  plane  angle  if  data  were  available  for 
confining  pressure  and  pore  water  pressure  characteristics  at  the  point  of  failure,  as  these 
variables  also  influence  failure  plane  angle. 

4.3.3  Critical  Geometry:  Darby-Thorne  Analysis 

In  section  4.3.2,  the  ability  of  each  of  the  four  tested  analyses  to  predict  incipient 
failure  block  geometry  was  assessed.  A  feature  of  computer  program  based  on  the  new 
analysis  developed  in  this  report  is  that  it  is  able  to  estimate  the  critical  geometry  of  initially 
stable  banks  by  applying  increments  of  lateral  fluvial  erosion  and  bed  degradation  to  the 
initial,  stable,  bank  profile,  until  critical  conditions  are  reached.  It  is,  therefore,  possible  to 
compare  predictions  of  failure  block  geometry  at  critical  conditions  with  the  observed  data 
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for  the  Darby-Thome  analysis.  This  is  a  useful  exercise,  since  the  failure  block  geometry  at 
the  point  of  failure  determines  the  amounts  of  bank  retreat,  channel  widening  and  volumes 
of  sediment  supplied  to  the  channel  once  the  banks  become  unstable.  To  determine  the 
ability  of  the  analysis  to  predict  these  parameters  at  critical  conditions  would  provide  insight 
into  the  utility  of  the  analysis  with  respect  to  modelling  channel  width  adjustments  due  to 
mass  failure  and  bank  erosion  processes. 

Table  5  Summary  of  calculated  mean  discrepancy  ratio  and  absolute  deviation  of  the 
discrepancy  ratio  for  critical  failure  block  geometry  conditions  using  the  Darbv-Thome 
bank  stability  analysis 


Variable 

Me 

Ad 

Failure  Plane  Angle,  (3 

0.993 

0.355 

Failure  Block  Volume,  VB 

2.703 

1.743 

Failure  Block  Width,  x 

1.430 

0.575 

Flood-plain  widening,  BW 

3.117 

1.681 

Critical  geometry  predictions  were  obtained  for  the  Darby-Thome  analysis  by 

running  the  analysis  using  observed  bank  geometries  and  "worst-case"  soil  properties.  If 

the  predicted  factor  of  safety  was  less  than  1.2,  the  failure  block  geometry  so  obtained  was 

assumed  to  be  the  predicted  critical  failure  block  geometry.  For  predicted  initial  factors  of 

safety  greater  than  1.2,  the  computer  program  based  on  this  analysis  was  used  to  supply 

increments  of  lateral  fluvial  erosion  and  bed  degradation  until  the  critical  conditions  were 

dW 

obtained.  For  each  analysis  a  ratio  of  0.1  was  assumed.  Predicted  critical  geometries  so 

obtained  were  then  compared  with  the  observed  data  (Figure  16  and  Table  5). 

It  can  be  seen  that  the  critical  failure  block  geometry  parameters  are  predicted  with 
varying  degrees  of  accuracy.  For  failure  plane  angle,  the  calculated  mean  discrepancy  ratio 
of  0.993  suggests  relatively  good  predictive  ability,  but  this  impression  is  somewhat 
misleading.  In  fact,  as  Figure  16a  indicates,  there  is  considerable  scatter  about  the  "line  of 
perfect  agreement".  When  compared  with  the  results  obtained  for  the  incipient  geometry 
analysis,  it  seems  that  low  failure  plane  angles  tend  to  be  raised  during  fluvial  erosion, 
while  the  steeper  failure  plane  angles  tend  to  remain  relatively  constant,  resulting  in 
improved  predictive  ability  overall,  but  with  much  increased  levels  of  uncertainty  for  any 
one  prediction.  This  tendency  is  as  expected,  since  the  steep  failure  planes  probably 
correspond  to  critical  banks,  so  that  these  banks  would  not  be  modified  during  the  critical 
geometry  analysis.  Alternatively,  fluvial  erosion  should  act  to  heighten  and  steepen  the 
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bank  profile  on  the  more  stable  banks,  which  in  turn  will  steepen  the  predicted  failure  plane 
angle. 

Failure  block  volumes  are  clearly  significantly  over-predicted,  with  a  calculated 
mean  discrepancy  ratio  of  2.703,  though  again  the  calculated  absolute  deviation  of  the 
discrepancy  ratio  and  Figure  16b  indicate  that  there  is  also  a  large  amount  of  variation  about 
the  "line  of  perfect  agreement".  The  increase  in  over-prediction  of  failure  block  volumes, 
compared  to  the  incipient  geometry  analysis,  corresponds  to  the  increasing  height  of  the 
failure  block  as  degradation  lowers  the  elevation  of  the  toe  of  the  bank. 

There  appears  to  be  little  difference  in  the  failure  block  width  parameters,  x  and 
BW,  predicted  by  the  incipient  and  critical  geometry  analyses,  indicating  that  the  shape  of 
the  critical  block  remains  relatively  constant  during  fluvial  erosion.  Thus,  increasing  failure 
plane  angles  tend  to  compensate  for  the  increases  in  bank  height,  to  maintain  relatively 
constant  block  widths  during  fluvial  erosion.  It  can  be  seen  that  total  block  width  has  a 
calculated  mean  discrepancy  ratio  of  1.430,  compared  to  the  flood-plain  widening 
increment  discrepancy  ratio  of  3.117.  Both  parameters  indicate  that  the  block  width  is 
significantly  over-predicted  by  the  Darby-Thome  analysis  at  critical  conditions. 

Overall,  it  appears  that  predictions  of  critical  geometry  obtained  using  the  Darby- 
Thome  analysis  do  not  compare  well  to  the  observed  data.  In  particular,  predictions  of 
block  magnitude  (volume  and  width)  are  significantly  over-predicted.  This  is  an  interesting 
result  with  regard  to  the  potential  for  applying  the  new  bank  stability  algorithm  in  the 
Darby-Thome  numerical  channel  evolution  model.  In  that  model,  under-prediction  of 
widening  rates  were  attributed  to  under-prediction  of  failure  block  width  and  neglect  of 
secondary  failures  (Darby  et  al..  In  Review).  Substituting  the  new  bank  stability  analysis 
into  the  channel  evolution  model  may  well  rectify  this  deficiency,  but  for  the  wrong 
reasons!  It  is  clear  that  estimates  of  land  loss  and  bank  sediment  yield  obtained  using  this 
new  analysis  may  be  significantly  in  error. 

On  the  other  hand,  these  data  in  isolation  can  be  somewhat  misleading.  It  must  be 
remembered  that  the  Darby-Thome  analysis  also  over-predicts  factor  of  safety,  and  that  the 
over-prediction  of  critical  failure  block  geometry  parameters  reflect  the  over-prediction  of 
factor  of  safety.  Moreover,  estimates  of  land  loss  and  bank  sediment  yield  obtained  using 
the  Darby-Thome  analysis  are  a  product  not  only  of  the  magnitude  of  individual  failures, 
but  also  the  frequency  at  which  those  failures  occur.  The  results  presented  in  this  report 
indicate  that  while  the  magnitude  of  block  size  will  be  over-predicted,  the  frequency  of 
failures  will  be  under-predicted  by  this  analysis.  The  result  of  this  counter-balancing  effect 
may  be  that  application  of  the  Darby-Thome  bank  stability  analysis  in  analyses  of  long-term 
morphological  response  may  well  result  in  successful  predictions,  though  it  is  recognised 
that  this  success  would  be  for  the  wrong  reasons.  The  important  point  is  that  it  is  clear  that 
the  Darby-Thome  analysis  is  the  most  physically-based,  most  widely  applicable  and  by  far 


52 


Predicted  Beta  (Degrees) 


most  successful  in  terms  of  predictive  ability  of  the  analyses  tested  here.  For  this  reason,  it 
is  recommended  that  the  bank  stability  analysis  should  be  adopted  in  studies  of  long  term 
river  channel  adjustments,  such  as  the  numerical  channel  evolution  model  developed  by 
Darby  and  Thome  (In  Review).  This  is  now  a  topic  of  current  research. 


5.  SUMMARY  AND  CONCLUSIONS 


The  ability  to  predict  the  stability  and  failure  geometry  of  eroding  riverbanks  is  an 
important  pre-requisite  for  modelling  numerically  alluvial  channel  adjustments,  as  well  as  a 
requirement  for  predicting  bank  erosion  rates  and  sediment  yield  associated  with  bank 
erosion.  However,  there  are  a  number  of  limitations  of  existing  bank  stability  analyses 
which  limit  their  physical  basis  and  predictive  ability.  Some  of  these  limitations  have  been 
addressed  here  through  the  development  of  a  new  bank  stability  analysis.  The  new 
approach  is  applicable  to  steep,  cohesive,  non-layered  riverbanks  which  fail  along  planar 
failure  surfaces.  The  analysis  is  developed  for  banks  with  profiles  characteristic  of  natural, 
eroding  riverbanks.  Failure  plane  angles  are  calculated  using  a  physically-based  approach; 
additionally  the  failure  plane  is  not  constrained  to  pass  through  the  toe  of  the  bank.  Pore 
water  and  hydrostatic  confining  pressure  effects  are  included  in  the  analysis.  Using  the 
procedure  developed  by  Darby  (1994),  the  new  analysis  also  is  capable  of  predicting  the 
probability  of  mass  failure,  as  well  as  an  overall  factor  of  safety  with  respect  to  mass 
failure. 

A  new  data  set,  comprising  information  about  51  individual  mass-failures,  has  been 
used  to  compare  the  predictive  abilities  of  4  bank  stability  analyses:  Lohnes  and  Handy 
(1968),  Huang  (1983),  Osman  and  Thorne  (1988)  and  the  analysis  developed  in  this 
report  Results  indicate  that  the  new  analysis  is  the  most  successful  of  the  tested  analyses  in 
terms  of  predicting  the  stability  of  riverbanks  with  respect  to  mass  failure.  However,  the 
new  analysis  still  significantly  over-predicts  bank  stabiUty.  Nevertheless,  it  is  demonstrated 
that  the  new  analysis  provides  significant  improvements  to  predictions  of  bank  stability 
when  compared  to  the  other  analyses.  It  is,  therefore,  concluded  that  application  of  the  new 
analysis  in  the  Darby-Thome  numerical  channel  evolution  model  (Darby,  1994;  Darby  and 
Thome,  In  Review;  Darby  et  al.,  In  Review)  would  result  in  improved  quantitative 
predictions  of  channel  adjustments.  Implementation  of  the  new  bank  stability  analysis  in  the 
numerical  channel  morphology  model  is  a  topic  of  current  research.  It  is  recommended  that 
the  new  analysis  should  also  be  used  as  a  stand  alone  tool  for  analysis  of  bank  stability  and 
bank  erosion  problems. 

The  data  set  compiled  to  test  the  various  bank  stability  analyses  is  also  a  resource  of 
considerable  value.  The  data  set  is  supplied  in  spreadsheet  format  on  a  diskette  with  this 
report.  The  data  set  is  the  first  suitable  for  testing  the  ability  of  bank  stability  analyses  to 
predict  accurately  the  factor  of  safety  and  geometry  of  individual  observed  mass-failures. 

The  new  bank  stability  analysis  may  be  used  as  a  stand-alone  analysis  to  aid 
engineers  in  the  design  of  stable  channels.  To  this  end,  the  analysis  has  been  developed  as 
a  computer  program,  which  is  supplied  on  diskette  with  this  report,  together  with 
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supporting  user  documentation.  The  computer  program  may  be  used  to  calculate  the  factor 
of  safety  or  probability  of  failure  of  riverbanks  of  geometry  specified  by  the  user.  The  user 
must  specify  the  geotechnical  characteristics  of  the  bank  material  and  select  whether  or  not 
to  include  pore  water  and/or  confining  pressure  effects  in  the  analysis.  If  the  user  does 
select  either  these  options  he  or  she  must  specify  the  ground  water  surface  elevation  and 
elevation  of  the  surface  water  in  the  channel,  respectively. 

A  particular  feature  of  the  computer  program  is  that  it  may  be  used  to  determine  the 
amounts  of  fluvial  erosion  and  bed  degradation  required  to  destabilize  an  initially  stable 
bank  to  the  point  of  failure.  This  feature  allows  the  computer  program  to  be  utilised  as  a 
design  aid  in  a  variety  of  stable  channel  design  applications,  such  as  determining  minimum 
grade-control  structure  elevations  required  to  prevent  degradation  initiating  channel 
widening  and  instability,  or  in  estimating  volumes  of  sediment  supplied  to  the  channel  from 
fluvial  erosion/bed  degradation  induced  bank  failures.  It  is  concluded  that  the  new  analysis 
and  computer  program  may  be  of  use  in  aiding  practicing  engineers  in  a  variety  of  problems 
which  require  analysis  of  bank  erosion  and  width  adjustments. 
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APPENDIX  II.  List  of  Symbols 


Ad  =  Absolute  deviation  of  the  discrepancy  ratio 
BW  =  Flood  plain  widening  increment 
BW*  =  "Potential"  flood  plain  widening  increment 
c  =  Soil  cohesion 

Fj.p  =  Hydrostatic  confining  pressure 
Fn  =  Function  in  equation  (39) 

Fn'  =  First  derivative  of  function  Fjj 
Fx  =  Horizontal  component  of  Fq, 

^xl  =  Horizontal  component  of  F^p  on  lowest  bank  face 
Fyl  =  Horizontal  component  of  F^p  on  middle  bank  face 
Fx3  =  Horizontal  component  of  F^p  on  upper  bank  face 
Fy  =  Vertical  component  of  Fq, 

FD  =  Driving  force  acting  on  incipient  failure  block 
FR  =  Resisting  force  acting  on  incipient  failure  block 
FS  =  factor  of  safety 
g  =  gravitational  acceleration 
H  =  Total  bank  height 
If  =  Uneroded  bank  height 

=  Effective  head  of  ground  water 
i  =  Bank  angle 
K  -  Tension  crack  depth 
/sTjj  =  Relic  tension  crack  depth 
L  =  Length  of  failure  surface 
Me  =  Mean  of  the  discrepancy  ratio 
n  =  Number  of  items  in  data  set 
AT  =  Computational  time  step 

=  Total  pore  pressure  acting  on  failure  plane 
=  Pore  pressure  acting  at  a  point  on  the  failure  plane 
V  =  Volume  of  failure  block  per  unit  reach  length 
Vft  =  Total  volumetric  bank  material  inflow  flux  due  to  mass  failure 

=  Reach-averaged  total  volumetric  bank  material  inflow  flux  due  to  mass  failure 
W  =  Weight  of  water  above  bank  face 
Vkf  =  Weight  of  incipient  failure  block 
X  =  Function  in  equation  (37) 

X'  =  First  derivative  of  X 
Xq  =  Observed  parameter 
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Xp  =  Predicted  parameter 
X  =  Total  failure  block  width 
X*  =  Potential  total  failure  block  width 
Ax  =  Model  reach  length 
Y  =  Function  in  equation  (38) 

T  =  First  derivative  of  Y 

3’f  =  Elevation  of  base  of  failure  plane 

y{*  =  Dimensionless  elevation  of  base  of  failure  plane 

jfp  =  Elevation  of  flood  plain 

=  Elevation  of  base  of  tension  crack 
3’s  =  Elevation  of  base  of  uneroded  bank  slope 
=  Elevation  of  base  of  relic  tension  crack 
yyf  =  Surface  water  elevation 

a  =  Angle  between  resultant  of  confining  pressure  and  normal  to  failure  plane 
P  =  Failure  plane  angle 
<j)  =  Soil  friction  angle 
Y=  Soil  unit  weight 
Yw  =  Unit  weight  of  water 
=  Density  of  water 

CO  =  Direction  of  resultant  of  confining  pressure 
q  =  Effective  head  of  surface  water 
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APPENDIX  III.  Bank  Geometry  Database  Users  Manual 


The  bank  geometry  database  used  to  test  the  predictive  abilities  of  the  selected  bank 
stability  analyses  in  this  study  is  supplied  on  two  3.5"  diskettes,  labelled  "Bank  Geometry 
Database  (Apple  Format)"  and  "Bank  Geometry  Database  (IBM  Format)",  respectively. 
The  diskettes  are  formatted  for  Apple  Macintosh  and  IBM  compatible  personal  computers, 
as  labelled.  For  the  Apple  Macintosh  disk,  data  is  presented  in  the  form  a  Microsoft  Excel 
(Version  3.0)  spreadsheet  package.  The  filename  for  the  data  base  is  "Geometry".  For  the 
IBM  compatible  disk,  data  is  presented  in  the  form  a  Lotus- 1,2, 3  (Version  4  for  Windows) 
spreadsheet  package.  The  filename  for  the  data  base  is  "Geometry". 

It  is  strongly  recommended  that  users  make  a  back-up  copy  of  the  files  containing 
the  data  base. 

The  database  consists  of  bank  geometry  and  soil  property  data  for  51  separate  mass 
failures  located  at  sites  along  eroding  creeks  in  northern  Mississippi.  The  data  is  supplied 
on  a  spreadsheet  so  that  users  can  easily  utilise  the  data  for  their  own  purposes.  Soil 
properties  for  "average"  and  "worst-case"  conditions  are  also  supplied.  Data  for  each  site  is 
arranged  in  rows  6  to  56,  inclusive,  in  the  spreadsheet,  with  individual  parameters 
containing  the  data  items  for  each  of  these  sites  arranged  in  separate  columns.  These 
columns  are  arranged  in  the  spreadsheet  in  the  following  way: 

Column  A.  Location  of  Site.  Text  data  in  this  column  provides  the  name  of  the  creek 
on  which  the  site  is  located. 

Column  B.  Data  Source.  Data  in  this  column  provides  a  reference  to  one  of  a  number 
of  footnotes  at  the  base  of  the  spreadsheet.  These  footnotes  provide  brief  details  concerning 
the  source  of  the  data  for  a  site  in  the  corresponding  row. 

Column  C.  Total  Bank  Height,  H  (m).  Provides  the  value  of  the  total  bank  height  at 
observed  critical  conditions. 

Column  D.  Uneroded  Bank  Height,  H'  (m).  Provides  the  value  of  the  uneroded 
bank  height  at  observed  critical  conditions. 

Column  E.  Bank  Angle  (Degrees).  Provides  the  value  of  the  uneroded  bank  slope  at 
observed  critical  conditions. 

Column  F.  Provides  the  tangent  of  the  bank  angle. 

Column  G.  Provides  the  sine  of  the  bank  angle. 

Column  H.  Provides  the  cosine  of  the  bank  angle. 

Column  I.  Tension  Crack  Depth,  K  (m).  Provides  the  value  of  the  tension  crack 
depth  observed  at  critical  conditions. 

Column  J.  Relic  Tension  Crack  Depth,  K^^  (m).  Provides  the  value  of  the  relic 
tension  crack  depth  observed  at  critical  conditions. 
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Column  K.  Failure  Plane  Angle,  P  (degrees).  Provides  the  value  of  the  failure 
surface  angle  observed  at  critical  conditions. 

Column  L.  Provides  the  tangent  of  the  failure  plane  angle. 

Column  M.  Provides  the  sine  of  the  failure  plane  angle. 

Column  N.  Provides  the  cosine  of  the  failure  plane  angle. 

Column  O.  Flood-Plain  Widening  Increment,  BW  (m).  Provides  the  value  of  the 
flood-plain  widening  increment  observed  at  the  point  of  failure. 

Column  P.  Total  Failure  Block  Width,  x  (m).  Provides  the  value  of  the  failure 
block  width  observed  at  the  point  of  failure. 

Column  Q.  Failure  Block  Volume,  VB  (m^m*^).  Provides  the  value  of  the  unit 
volume  of  the  failure  block  at  the  point  of  failure. 

Column  R.  Blank  Column. 

Column  S.  Average  Soil  Unit  Weight,  y  (kNm*^).  Provides  the  value  of  the 
average,  or  "best  case"  soil  unit  weight.  This  corresponds  to  the  measured  value  under 
"dry"  conditions. 

Column  T.  Average  Soil  Cohesion,  c  (kPa).  Provides  the  value  of  the  average,  or 
"best  case"  soil  cohesion.  This  corresponds  to  the  measured  value  under  "dry"  conditions. 
The  method  of  measurement  is  referenced  in  Column  W. 

Column  U.  Average  Soil  Friction  Angle,  ([)  (degrees).  Provides  the  value  of  the 
average,  or  "best  case"  soil  friction  angle.  This  corresponds  to  the  measured  value  under 
"dry"  conditions.  The  method  of  measurement  is  referenced  in  Column  W. 

Column  V.  Provides  the  tangent  of  soil  friction  angle. 

Column  W.  References  the  method  of  measurement  of  the  soil  properties.  The  characters 
in  this  column  are  cross-referenced  with  a  foot  note  at  the  base  of  the  spreadsheet. 

Column  X.  Worst  Case  Soil  Unit  Weight,  y  (kNm'^).  Provides  the  value  of  the 
"worst  case"  soil  unit  weight.  This  corresponds  to  the  value  estimated  for  "saturated" 
conditions. 

Column  Y.  Worst  Case  Soil  Cohesion,  C  (kPa).  Provides  the  value  of  the  "worst 
case"  soil  cohesion.  This  corresponds  to  the  value  estimated  for  "saturated"  conditions.  The 
method  of  measurement  is  referenced  in  Column  AB. 

Column  Z.  Worst  Case  Soil  Friction  Angle,  (j)  (degrees).  Provides  the  value  of 
the  "worst  case"  soil  friction  angle.  This  corresponds  to  the  value  estimated  for  "saturated" 
conditions.  The  method  of  measurement  is  referenced  in  Column  AB. 

Column  AA.  Provides  the  tangent  of  the  "worst-case"  soil  friction  angle. 

Column  AB.  References  the  method  of  measurement  of  the  soil  properties.  The 
characters  in  this  column  are  cross-referenced  with  a  foot  note  at  the  base  of  the 
spreadsheet. 

Column  AC.  Blank  Column. 
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Column  AD.  Total  Pore  Pressure,  (kPa).  This  data  is  unavailable. 

Column  AE.  Confining  Pressure,  F^p  (kPa).  This  data  is  unavailable. 

Column  AF.  Pore  Pressure  ratio,  ru.  This  data  is  estimated  for  "worst  case" 
conditions  for  saturated  soils  (see  Huang,  1983). 

Column  AG.  Location  of  Failure  Plane,  yf*.  This  gives  the  dimensionless  bank 
height  (measured  in  an  upwards  direction  from  the  base  of  the  bank)  at  which  the  most 
critical  failure  surface  is  located. 
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APPENDIX  IV.  Bank  Stability  Model  Users  Manual 


Requirements  to  Use  the  Computer  Program 

The  source  code  of  the  computer  model  based  on  the  bank  stability  analysis 
developed  in  this  report  is  written  in  FORTRAN-77  and  is  contained  on  an  MS-DOS 
formatted  disk,  labelled  "Bank  Program",  along  with  an  executable  version  of  the  code. 
This  disk  is  supplied  with  original  copies  of  this  report  The  disk  also  holds  4  example  data 
files  required  to  run  the  model:  UNITWEIG.DAT,  COHESION.DAT,  FRICTION.DAT 
and  SOILDATA.DAT.  The  source  code  is  titled  DARBY.FOR  and  the  executable  code  is 
titled  DARBY.EXE.  To  run  the  computer  program,  an  IBM  compatible  PC  with  math  co¬ 
processor  installed  is  required.  It  is  recommended  that  users  should  make  back-up  copies 
of  the  6  files  contained  on  the  diskette,  and  load  these  files  onto  their  PCs.  It  is 
recommended  that  the  computer  program  and  supporting  data  files  should  be  installed  in  a 
single  directory  on  the  C:  drive  of  the  PC.  Users  may  freely  use  the  program  supplied  with 
this  report,  but  may  not  distribute  the  code  to  other  parties.  The  authors  accept  no 
responsibility  or  liability  resulting  from  use  of  the  computer  program. 

About  the  Manual 

In  this  Appendix,  details  are  provided  which  explain  how  to  enter  data  to  the  model, 
how  to  select  various  user-specified  options,  and  how  to  read  the  data  stored  in  the  output 
file.  Worked  examples,  using  the  supplied  example  data  files,  are  also  presented  to 
demonstrate  the  application  of  the  computer  program.  In  using  this  manual,  it  is  assumed 
that  the  user  has  read  and  understood  the  main  body  of  this  report  Together,  the  report  and 
manual  provide  sufficient  information  for  users  to  train  themselves  to  operate  the  computer 
program  and  apply  it  to  analyse  bank  stability  problems. 

A.  ENTERING  DATA  TO  THE  MODEL 

The  bank  stability  analysis  developed  in  this  report  requires  data  relating  to  the 
geometry  of  the  bank  profile,  the  geotechnical  properties  of  the  bank  soil,  as  well  as 
information  relating  to  the  ground  and  surface  water  elevations.  This  data  is  supplied  to  the 
computer  program  in  two  ways.  First,  some  data  is  stored  in  input  data  files,  which  must 
be  edited  prior  to  execution  of  the  program.  Second,  data  is  supplied  to  the  model  directly, 
in  response  to  a  range  of  prompts  that  appear  on  the  user's  screen  during  model  execution. 
It  should  be  noted  that  all  data  is  entered  to  the  model  in  SI  units,  as  described  in  the  main 
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report.  Reference  should  be  made  to  the  report  if  there  is  any  doubt  about  the  appropriate 
units  to  use  for  a  particular  input  parameter. 

1.  Editing  Input  Data  Files 

It  is  necessary  to  check  and  edit  a  total  of  4  input  data  files  prior  to  running  the 
computer  program.  The  file  UNITWEIG.DAT  contains  information  about  the  frequency  of 
occurrence  of  various  values  of  soil  unit  weight  for  the  analysed  bank  material.  Similarly, 
the  files  FRICTION.DAT  and  COHESION.DAT  contain  information  relating  to  the 
frequency  of  occurrence  of  the  values  of  soil  friction  angle  and  cohesion,  respectively,  for 
the  bank  material.  These  three  files  are  used  to  supply  data  required  for  the  probabilistic 
bank  stability  analysis.  The  file  SOILDATA.DAT  contains  summary  information  relating  to 
the  soil  property  frequency  distributions  contained  in  the  previous  three  files.  This  file  is 
used  primarily,  but  not  exclusively,  to  supply  data  required  for  the  non-probabilistic  bank 
stability  analysis. 

Between  them  these  files  contain  all  the  information  relating  to  the  geotechnical 
properties  of  the  bank  soils  necessary  to  execute  the  computer  program.  Users  are  warned 
that  great  care  must  be  taken  in  editing  the  input  data  files.  Care  must  be  taken  to  ensure  that 
not  only  are  correct  values  entered  for  individual  parameters,  but  that  the  data  in  each  file  is 
entered  in  the  correct  format.  Also,  certain  constraints,  in  addition  to  those  outlined  in 
section  3.9  of  the  main  report,  exist  on  the  range  of  values  of  some  variables.  These  issues 
are  now  addressed  in  turn  for  each  of  the  input  data  files. 

UNITWEIG.DAT 

This  file  contains  data  regarding  the  probability  of  occurrence  of  particular  values  of  the  soil 
unit  weight.  Data  are  entered  in  two  columns.  An  example  of  the  format  of  this  file  is 
provided  in  Table  AIVl,  which  corresponds  to  the  file  UNITWEIG.DAT  supplied  on  the 
diskette.  The  first  column  of  numbers  in  the  file  represents  the  soil  unit  weight  data,  in  the 
form  of  positive,  non-zero  values.  Note  that  the  soil  unit  weight  values  are  entered  in  Nm' 
,  not  the  more  commonly  used  kNm'-^.  The  second  column  of  numbers  corresponds  to  the 
probability  of  occurrence  (non-dimensional  units)  in  the  range  (0  <  P(y)  <  1)  of  each  of  the 
unit  weight  values.  Since  the  second  column  denotes  probabilities  of  occurrence,  the  values 
entered  in  this  column  must  sum  to  a  total  of  unity.  Note  that  the  first  field  is  left  blank  (x 
denotes  a  blank  field).  Soil  unit  weight  data  entered  in  the  first  column  then  must  occupy 
the  next  7  fields,  which  includes  one  decimal  place.  Two  fields  are  left  blank  to  separate  the 
first  and  second  columns.  The  probability  data  in  the  second  column  must  occupy  the  next 
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5  fields,  including  3  decimal  plates.  Note  that  numeric  characters  must  be  entered  into  every 
field  of  both  columns  (e.g.  07000.0  must  be  entered,  not  7000  or  7000.0). 


Table  AT VI  Example  of  format  of  UNITWEIG.DAT  input  data  file 


Field 

1 

2 

3 

m 

5 

6 

m 

8 

9 

0 

1 

2 

3 

m 

5 

Soil  Unit  Weight  Value  (Nm‘3) 

Blank 

Probability 

X 

0 

Di 

0 

0 

0 

■ 

0 

X 

X 

0 

■ 

0 

0 

0 

X 

0 

9 

0 

0 

0 

■ 

0 

X 

X 

0 

■ 

0 

0 

0 

X 

1 

1 

0 

0 

0 

0 

X 

X 

0 

■ 

0 

0 

0 

X 

1 

3 

0 

0 

0 

HHIII 

0 

X 

X 

0 

1 

2 

5 

X 

1 

5 

0 

0 

0 

■ 

0 

X 

X 

0 

■ 

0 

6 

2 

1 

m 

0 

0 

0 

0 

X 

X 

0 

■ 

1 

8 

8 

91 

1 

9 

0 

0 

0 

■ 

0 

X 

X 

0 

6 

2 

5 

X 

2 

1 

0 

0 

0 

0 

X 

X 

0 

0 

0 

0 

X 

2 

3 

0 

0 

0 

■ 

0 

X 

X 

0 

. 

0 

0 

0 

Note  that  the  total  number  of  rows  in  the  file  is  designated  as  a  positive  integer  called 
NUM2.  This  integer  is  later  required  as  input  data  in  the  file  SOILDATA.DAT.  In  the 
above  example  NUM2  =  9.  The  maximum  number  of  rows  permitted  is  59,  so  that  the 
condition  1  <  NUM2  <  59  must  be  satisfied.  The  file  UNITWEIG.DAT  is  not  used  in  the 
non-probabilistic  analysis.  If  it  is  intended  to  use  only  the  non-probabilistic  analysis,  it  is, 
therefore,  not  necessary  to  edit  this  file. 

FRICTION.DAT 

This  file  contains  data  regarding  the  probability  of  occurrence  of  particular  values  of  the  soil 
friction  angle.  Data  are  again  entered  in  two  columns.  An  example  of  the  format  of  this  file 
is  illustrated  in  Table  AIV2,  which  is  the  file  FRICTION.DAT  supplied  on  the  disk  with 
this  report.  The  first  column  of  numbers  in  the  file  represents  soil  friction  angle  data,  in  the 
form  of  positive,  non-zero  values.  Soil  friction  angle  values  are  entered  in  degrees.  The 
second  column  of  numbers  corresponds  to  the  probability  of  occurrence  (non-dimensional 
units)  in  the  range  (0  <  P((())  <  1)  of  the  corresponding  friction  angle  value.  Since  the 
second  column  denotes  probabilities  of  occurrence,  the  values  entered  in  this  column 
should  always  sum  to  a  total  of  unity.  Note  that  the  first  field  is  left  blank  (x  denotes  a 
blank  field).  Soil  friction  angle  data  entered  in  the  first  column  then  must  occupy  the  next  4 
fields,  which  includes  one  decimal  place.  Two  fields  are  left  blank  to  separate  the  first  and 
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second  columns.  The  probability  data  in  the  second  column  must  occupy  the  next  5  fields, 
including  3  decimal  places.  Note  that  numerical  characters  must  be  entered  into  every  field 
of  both  columns  (e.g.  07.5  is  entered,  not  7.5). 


Table  AIV2  Example  of  format  ofFRICTION.DAT  input  data  file 


Field 

1 

2 

3 

m 

5 

6 

7 

8 

9 

0 

1 

2 

Friction  Angle  (degs) 

Blank 

Probability 

X 

0 

m 

5 

X 

X 

0 

0 

0 

0 

X 

1 

2 

. 

5 

X 

X 

0 

0 

6 

3 

X 

1 

m 

■ 

5 

X 

X 

0 

■ 

0 

0 

0 

X 

2 

2 

■ 

5 

X 

X 

0 

■ 

1 

2 

5 

X 

2 

m 

5 

X 

X 

0 

1 

2 

5 

X 

3 

2 

. 

5 

X 

X 

0 

3 

m 

5 

X 

3 

m 

■ 

5 

X 

X 

0 

2 

5 

0 

X 

■ 

2 

5 

X 

X 

0 

• 

0 

6 

2 

Note  that  the  total  number  of  rows  in  the  file  is  designated  as  a  positive  integer  called 
NUM3.  This  integer  is  later  required  as  input  data  in  the  file  SOILDATA.DAT.  In  the 
above  example  NUM3  =  8.  The  maximum  number  of  rows  permitted  is  59  (i.e.  1  <  NUM3 
<  59).  The  file  UNITWEIG.DAT  is  not  used  in  the  non-probabilistic  analysis.  If  it  is 
intended  to  use  only  the  non-probabilistic  analysis,  it  is,  therefore,  not  necessary  to  edit  this 
file. 

COHESION.DAT 

This  file  contains  data  regarding  the  probability  of  occurrence  of  particular  values  of  the  soil 
cohesion.  Data  are  entered  in  two  columns.  An  example  of  the  format  of  this  file  is 
provided  in  Table  ArV3,  which  corresponds  to  the  file  COHESION.DAT  supplied  on  the 
diskette.  The  first  column  of  numbers  in  the  file  represents  the  soil  cohesion  data,  in  the 
form  of  positive,  non-zero  values.  Note  that  the  soil  cohesion  values  are  entered  in  Nm'^, 
not  the  more  commonly  used  kNm‘^.  The  second  column  of  numbers  corresponds  to  the 
probability  of  occurrence  (non-dimensional  units)  in  the  range  (0  <  P(y)  ^  1)  of  each  of  the 
cohesion  values.  Since  the  second  column  denotes  probabilities  of  occurrence,  the  values 
entered  in  this  column  must  sum  to  a  total  of  unity.  Note  that  the  first  field  is  left  blank  (x 
denotes  a  blank  field).  Soil  cohesion  data  entered  in  the  first  column  then  must  occupy  the 
next  7  fields,  which  includes  one  decimal  place.  Two  fields  are  left  blank  to  separate  the 
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first  and  second  columns.  The  probability  data  in  the  second  column  must  occupy  the  next 
5  fields,  including  3  decimal  plates.  Note  that  numeric  characters  must  be  entered  into  every 
field  of  both  columns  (e.g.  02500.0  must  be  entered,  not  2500  or  2500.0). 


Table  ATV3  Example  of  format  of  COHESION.DAT  input  data  file 


Field 

1 

2 

3 

m 

5 

6 

m 

8 

9 

0 

1 

2 

3 

m 

5 

Soil  Cohesion  Value  (Nm"^) 

Blank 

Probability 

X 

0 

2 

5 

0 

0 

0 

X 

X 

0 

0 

0 

0 

X 

0 

m 

5 

0 

0 

0 

X 

X 

0 

1 

8 

7 

X 

1 

2 

5 

0 

0 

. 

0 

X 

X 

0 

. 

2 

5 

0 

X 

1 

m 

5 

0 

0 

0 

X 

X 

0 

■ 

m 

3 

8 

X 

2 

2 

5 

0 

0 

■ 

0 

X 

X 

0 

0 

0 

0 

X 

2 

Di 

5 

0 

0 

. 

0 

X 

X 

0 

■ 

0 

0 

0 

X 

3 

2 

5 

0 

0 

■ 

0 

X 

X 

0 

1 

2 

5 

X 

3 

m 

5 

0 

0 

0 

X 

X 

0 

. 

0 

0 

0 

X 

■ 

2 

5 

0 

0 

■ 

0 

X 

X 

0 

■ 

0 

0 

0 

Note  that  the  total  number  of  rows  in  the  file  is  designated  as  a  positive  integer  called 
NUM4.  This  integer  is  later  required  as  input  data  in  the  file  SOILDATA.DAT.  In  the 
above  example  NUM4  =  9.  The  maximum  number  of  rows  permitted  is  59,  so  that  the 
condition  1  <  NUM4  <  59  must  be  satisfied.  The  file  COHESION.DAT  is  not  used  in  the 
non-probabilistic  analysis.  If  it  is  intended  to  use  only  the  non-probabilistic  analysis,  it  is, 
therefore,  not  necessary  to  edit  this  file. 

SOILDATA.DAT 

This  file  contains  the  modal  values  of  the  soil  cohesion,  friction  angle  and  unit  weight 
frequency  distributions,  which  are  required  to  run  the  non-probabilistic  bank  stability 
analysis.  If  the  frequency  distributions  are  unavailable  representative  values  of  these  soil 
parameters  may  be  substituted  for  the  modal  values.  In  addition  the  file  contains  the 
integers  NUM2,  NUM3  and  NUM4  which  are  used  in  the  probabilistic  bank  stability 
analysis,  and  which  refer  to  the  number  of  rows  entered  in  the  UNITWEIG.DAT, 
FRICTION.DAT  and  COHESION.DAT  data  files,  respectively.  It  is,  therefore,  necessary 
to  check  this  input  data  file  whichever  bank  stability  analysis  is  selected. 

Data  is  input  to  the  file  along  1  row  in  6  columns,  which  occupy  a  variety  of  fields. 
An  example  of  the  format  of  this  file  is  provided  in  Table  AIV4,  which  corresponds  to  the 
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file  SOILDATA.DAT  supplied  on  the  diskette.  The  first  field  is  left  blank.  The  first  column 
of  numeric  characters  represents  the  cohesion,  entered  in  Pa.  This  column  occupies  7 
fields,  of  which  1  must  be  a  decimal  place.  Two  blank  fields  separate  the  first  and  second 
columns.  Numeric  characters  in  the  second  column  represents  the  soil  unit  weight  in  Nm’^. 
The  second  column  also  occupies  7  fields,  of  which  1  must  be  a  decimal  place.  The  second 
and  third  columns  are  separated  by  two  blank  fields.  Numeric  characters  in  the  third 
column  represents  the  soil  friction  angle,  in  degrees.  This  column  occupies  4  fields,  1  of 
which  must  be  a  decimal  place.  Two  blank  fields  separate  the  third  and  fourth  columns. 
Numeric  characters  entered  in  the  fourth  column  represents  the  integer  NUM2.  This  value 
occupies  3  fields.  Finally,  the  integers  NUM3  and  NUM4  occupy  three  fields  each  in  the 
fifth  and  sixth  columns,  respectively.  Again,  these  columns  are  separated  by  two  blank 
fields. 

Table  AIV4  Example  of  format  of  SOILDATA.DAT  int?ut  data  file 

xl7500.0xxl9000.0xx32.5xx009xx008xx009 

WARNING 

When  entering  data  into  the  input  data  files,  users  are  warned  that  great  care  must  be  taken 
to  follow  the  procedures  outlined  above.  In  addition,  as  explained  in  the  main  report  (see 
section  3.9),  great  care  must  be  taken  to  specify  realistically  the  soil  property  frequency 
distributions  when  applying  the  probabilistic  bank  stability  analysis.  Particular  care  must  be 
paid  in  specifying  soil  unit  weight  distributions  which  have  low  magnitudes  of  soil  unit 
weight  with  non-zero  probabilities.  When  using  the  pore  water  pressure  analysis,  it  is 
possible  that  some  combinations  of  bank  geometry  and  pore  water  pressure  values,  coupled 
with  occurrence  of  low  magnitude  soil  unit  weight  values,  may  lead  to  instances  where 
computed  failure  block  weights  are  less  than  computed  pore  water  pressures.  This  is 
physically  meaningless  and  spurious  numerical  results  may,  therefore,  be  obtained  in  those 
circumstances. 

2.  Entering  Data  at  the  Terminal 

After  editing  the  input  data  files,  program  execution  commence  by  typing  the  name 
of  the  executable  code,  DARBY,  and  hitting  the  RETURN  key.  A  sequence  of  prompts 
will  appear  requesting  the  user  to  input  various  data  items.  These  data  items  can  be  sub¬ 
divided  into  a  number  of  categories.  These  are  now  described  in  the  order  in  which  they 
appear  on  the  screen. 
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2.1  Bank  Geometry  Data 


Data  prompts  in  this  category  request  the  various  data  which  describe  the  shape  of 
the  bank  profile  (see  the  definition  diagram  Figure  4  in  the  main  report).  Users  should  note 
that  all  elevation  data  must  be  entered  in  metres,  and  should  be  input  in  absolute  elevations 
(e.g.  metres  above  sea  level  or  other  datum).  In  sequence,  the  following  prompts  appear: 

"Enter  flood-plain  elevation  (yfp)" 

At  this  prompt,  enter  the  elevation  of  the  flood-plain  (in  metres),  as  defined  by  the  point  yfp 
in  Figure  4.  Then  hit  the  RETURN  key  to  obtain  the  next  prompt: 

"Enter  elevation  of  bank  toe  (yf)" 

Type  the  elevation  of  the  bank  toe  (in  metres),  as  defined  by  the  point  yf  in  Figure  4.  Then 
hit  the  RETURN  key  to  obtain  the  next  prompt: 

"Enter  bank  angle  (degrees)" 

Type  the  bank  angle  (in  degrees),  defined  by  the  angle  i  in  Figure  4.  The  value  of  the  bank 
angle  should  be  in  the  range  50  <  i  <  90  to  ensure  computational  integrity  (see  section  3.9 
of  main  report).  Then  hit  the  RETURN  key  to  obtain  the  next  prompt: 

"Enter  elevation  of  base  of  bank  slope  (ys)" 

Type  the  elevation  of  the  base  of  the  uneroded  bank  slope  (in  metres)  defined  by  the  point 
y^  in  Figure  4.  Note  that  if  the  bank  is  completely  uneroded  then  y^  =  yf  should  be  used. 
Then  hit  the  RETURN  key  to  obtain  the  next  prompt: 

"Enter  elevation  of  base  of  vertical  slope  (yt)" 

Type  the  elevation  (in  metres)  of  the  base  of  any  relic  tension  crack  that  may  be  present, 
defined  by  the  point  yj  in  Figure  4.  If  there  is  no  tension  crack  present  use  y^  =  yfp.  Then 
hit  the  RETURN  key  to  obtain  the  next  prompt: 
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•’Enter  elevation  of  base  of  tension  crack  (yk)" 


Type  the  elevation  (in  metres)  of  the  base  of  the  tension  crack  (if  present),  defined  by  the 
point  y|j  in  Figure  4.  Users  are  reminded  that  values  of  tension  crack  depth  greater  than  half 
the  total  bank  height  may  lead  to  negative  values  of  failure  block  weight,  which  is 
impossible.  The  computer  program  automatically  constrains  tension  cracks  so  they  do  not 
exceed  this  maximum  depth.  This  should  be  considered  when  interpreting  analyses  for 
which  deeper  cracks  were  (incorrectly)  specified.  If  there  is  no  tension  crack  present  use 
=  yfp.  Then  hit  the  RETURN  key  to  obtain  the  next  prompt,  which  is  in  the  next  category. 

2.2  User-Specified  Output  File  Name 
The  prompt: 

"Please  name  the  output  data  file" 

requests  the  user  to  name  the  output  data  file  to  which  the  results  of  the  analysis  will  be 
written.  The  user  may  use  any  combination  of  letters  or  numeric  characters,  up  to  a 
maximum  of  10  characters.  In  addition  to  these  10  characters,  the  named  output  data  file 

will  automatically  be  given  the  extension _ .DAT.  For  example,  if  the  user  names 

the  file  "TESTl",  the  complete  output  file  name  will  be  named  "TEST1.DAT".  Users 
should  not  attempt  to  add  this  extension  to  their  file  name.  After  typing  the  file  name, 
hitting  the  RETURN  key  obtains  the  next  prompt. 

2.3  Critical  Geometry  Analysis 
The  prompt 

"Do  you  want  to  find  the  critical  geometry  (y/n)" 

allows  the  user  to  select  whether  or  not  to  analyse  only  the  stability  of  the  initially  specified 
bank  profile,  or  whether  to  use  the  program  to  estimate  the  amount  of  fluvial  erosion 
required  to  bring  that  initially  specified  bank  to  the  point  of  failure.  If  the  critical  geometry 
analysis  is  required,  the  user  should  type  a  lower-case  "y"  and  hit  the  RETURN  key,  to 
obtain  the  next  prompt.  Typing  any  other  key,  including  an  upper-case  "Y"  will  result  in 
the  critical  geometry  analysis  not  being  selected.  It  is  stressed  that  only  a  lower  case  "y"  is 
permitted  if  the  user  wishes  to  select  this  analysis.  It  should  also  be  noted  that  if  this 
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analysis  is  selected,  the  user  must  also  select  the  non-probabilistic  bank  stability  analysis 
(see  below).  If  the  analysis  is  selected  the  next  prompt: 

"Enter  bank  erosion  parameter  (dW/dZ)" 

allows  the  user  to  specify  the  ratio  of  lateral  erosion  to  bed  degradation  supplied  to  the  bank 
profile.  In  the  absence  of  field  data,  it  is  suggested  that  the  value  of  0.1  may  be  entered  as  a 
reasonable  default  value  for  many  cohesive  riverbanks  not  subjected  to  severe  particle-by- 
particle  erosion.  Hitting  the  RETURN  key  obtains  the  next  prompt. 

2.4  Ground  Water  Data 

The  next  prompt  allows  the  user  to  select  whether  or  not  to  include  pore  water 
pressure  effects  in  the  bank  stability  analysis.  In  response  to  the  prompt: 

"Do  you  want  to  include  pore  pressure  (y/n)" 

the  user  should  type  a  lower-case  "y"  and  hit  RETURN  only  if  pore  water  pressure  effects 
are  to  be  included  in  the  bank  stability  analysis.  Typing  any  other  key  will  result  in  the  pore 
water  pressure  analysis  not  being  selected.  If  the  analysis  is  selected,  the  prompt: 

"Enter  ground- water  surface  elevation  (gwse)" 

requests  the  elevation  (in  metres)  of  the  ground  water  surface,  defined  by  the  line  GWSE  in 
Figure  4.  Any  positive  value  in  the  range  (0  <  GWSE  <  y^)  may  be  entered,  though  recall 
that  analyses  with  values  close  to  the  flood-plain  elevation  may  encounter  difficulties  for 
certain  bank  geometries  and  low  values  of  soil  unit  weight  (probabilistic  analysis).  Then  hit 
the  RETURN  key  to  obtain  the  next  prompt. 

2.5  Surface  Water  Data 

The  next  prompt  allows  the  user  to  select  whether  or  not  to  include  hydrostatic 
confining  pressure  effects  in  the  bank  stability  analysis.  In  response  to  the  prompt: 

"Do  you  want  to  include  pore  pressure  (y/n)" 

the  user  should  type  a  lower-case  "y"  and  hit  RETURN  only  if  hydrostatic  confining 
pressure  effects  are  to  be  included  in  the  bank  stability  analysis.  Typing  any  other  key  will 
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result  in  the  confining  pressure  analysis  not  being  selected.  If  the  analysis  is  selected,  the 
prompt: 


"Enter  surface  water  elevation  (wse)" 

requests  the  elevation  (in  metres)  of  the  surface  water  surface,  defined  by  the  line  WSE  in 
Figure  4.  Any  positive  value  in  the  the  range  (0  <  WSE  <  y^p)  may  be  entered.  Then  hit  the 
RETURN  key  to  obtain  the  next  prompt. 

2.6  Selection  of  Bank  Stability  Analysis 
The  prompt: 

"Choose  modal  (np)  or  probabilistic  (p)  analysis" 

allows  the  user  to  select  between  either  of  the  non-probabilistic  (modal)  or  probabilistic 
bank  stability  analyses.  If  the  user  has  previously  selected  the  critical  geometry  analysis, 
the  user  must  select  the  non-probabilistic  analysis  at  this  point.  To  do  this,  the  user  should 
type  lower-case  "np"  and  hit  the  RETURN  key.  Typing  any  other  key  will  select  the 
probabilistic  bank  stability  analysis.  In  either  case,  hitting  RETURN  obtains  the  next 
prompt. 

Users  should  note  that  when  the  non-probabilistic  analysis  is  selected,  soil  input 
data  from  the  file  SOILDATA.DAT  is  used.  The  probabilistic  analysis  utilises  soil  data 
from  the  COHESION.DAT,  UNITWEIG.DAT  and  FRICTION.DAT  files.  The 
appropriate  files  should  be  checked  carefully  when  interpreting  the  output  data. 

2.7  Data  Check  Prompt 
The  final  prompt: 

"Please  check  data  and  options  (type  ok  to  proceed)" 

allows  the  user  to  check  and,  if  necessary,  correct  the  data  input  in  response  to  the 
prompts.  If  no  amendments  are  necessary,  type  lower-case  "ok"  and  hit  RETURN  to 
continue  program  execution.  If  the  user  wishes  to  amend  any  user-specified  data  items, 
type  any  other  key  and  hit  RETURN.  The  data  input  process  is  then  repeated  from  the  first 
prompt.  Note  that  file  input  data  cannot  be  amended  at  this  point.  If  an  error  is  spotted  in 
the  input  data  files,  it  will  be  necessary  to  terminate  program  execution. 
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B.  PROGRAM  EXECUTION 


After  selecting  the  proceed  option  in  response  to  the  data  check  prompt,  hitting  the 
RETURN  key  results  in  continued  program  execution.  Processing  will  continue  until 
either: 

1.  A  message  is  displayed  indicating  program  termination.  This  indicates  successful 
program  execution.  The  user  can  now  interpret  the  output  data  (see  below). 

2.  A  program-supplied  diagnostic  error  message  is  displayed,  followed  by  the  termination 
message.  The  program- supplied  diagnostic  message  will  indicate  the  possible  reason  for 
program  failure  (see  "error  messages",  below). 

3.  A  compiler-supplied  error  message  is  displayed,  followed  by  the  termination  message. 
Compiler-supplied  error  messages  do  not  usually  indicate  the  reason  for  failure,  other  than 
perhaps  locating  the  point  in  the  source  code  at  which  the  program  "crashed".  In  the 
unlikely  event  that  such  a  failure  occurs,  it  is  recommended  that  the  user  should  carefully 
check  all  input  data.  If  this  check  fails  to  locate  the  reason  for  the  error,  contact  the  authors 
directly  at  the  address  on  the  cover  of  this  report. 

C.  OUTPUT  DATA 

Once  program  execution  is  complete,  output  data  is  written  to  the  file _ .DAT 

previously  named  by  the  user.  To  look  at  the  output  data,  users  must  open  their  output  file 
using  an  editor.  In  the  file,  data  is  presented  in  two  sections;  an  input  data  summary  and  an 
output  data  summary.  The  input  data  summary  describes  the  soil  properties  entered  into  the 
model,  together  with  messages  outlining  the  various  options  selected  by  the  user.  The  soils 
data  presented  in  the  input  data  summary  correspond  to  the  modal  soil  values  entered  in  the 
file  SOILDATA.DAT.  These  values  are  used  in  the  non-probabilistic  bank  stability 
analysis,  but  also  serve  to  provide  a  summary  of  the  representative  soil  properties  specified 
in  the  probabilistic  bank  stability  analysis.  In  interpreting  the  results,  care  should  be  taken 
that  the  modal  soil  values  are  accurately  cross-referenced  with  the  soil  property  frequency 
distributions. 

The  output  data  summary  details  the  predicted  stability  of  the  specified  bank, 
together  with  information  about  the  geometry  of  the  most-critical  predicted  failure  surface. 
If  the  critical  geometry  analysis  was  selected,  the  degradation  required  to  cause  failure  and 
the  height  of  the  critical  bank  is  also  displayed.  Information  is  provided  relating  to  the 
factor  of  safety  (non-probabilistic  analysis)  or  probability  of  failure  (probabilistic  analysis) 
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of  the  most  critical  bank.  In  the  non-probabilistic  analysis,  the  bank  is  predicted  to  be 
unstable  with  respect  to  mass  failure  if  the  factor  of  safety  is  less  than  or  equal  to  unity. 
Bank  stability  increases  as  the  factor  of  safety  increases  above  unity,  or  as  the  probability 
of  failure  decreases.  Output  data  relating  to  details  of  the  geometry  of  the  most  critical 
failure  surface,  as  defined  in  Figure  4,  are  then  provided,  for  both  types  of  analysis.  In 
order,  the  output  data  summary  lists  the  most  critical  failure  plane  angle,  P  (degrees),  most 
critical  failure  block  volume,  VB  (m^/m),  total  failure  block  width,  x,  (m),  flood-plain 
widening  increment,  BW  (m)  and  dimensionless  height  of  the  most  critical  failure  surface 
(Vf*).  This  latter  parameter  is  non-dimensionalized  by  dividing  the  height  of  the  location  of 
the  most  critical  failure  surface  above  the  base  of  the  bank  by  the  total  bank  height.  Thus  a 
basal  failure  corresponds  to  yf*  =  0. 

D.  ERROR  MESSAGES 

A  number  of  program-supplied  diagnostic  error  messages  have  been  built  into  the 
program,  to  indicate  probable  reasons  for  program  failure.  These  messages  are  now  listed. 

1.  Warning  -  pore  water  pressure  options  failed.  This  message  indicates  that  the 
program  failed  to  select  any  of  the  possible  solutions  for  pore  water  pressure.  Check  that 
realistic  values  of  jfp,  yf,  ys»  )’t’  GWSE  and  WSE  have  been  selected. 

2.  Warning  -  confining  pressure  option  failed.  This  message  indicates  that  the 
program  failed  to  select  any  of  the  possible  solutions  for  hydrostatic  confining  pressure. 
Check  that  realistic  values  of  yfp,  yf,  yg,  y^,  y^  GWSE  and  WSE  have  been  selected. 

3.  Warning  -  failure  plane  options  failed.  This  message  indicates  that  the  program 
failed  to  select  any  of  the  possible  solutions  for  failure  plane  angle.  Check  that  realistic 
values  of  yfp,  yf,  y^,  y^.,  y^  GWSE  and  WSE  have  been  selected. 

4.  error  -  computed  probability  is  more  than  1.0.  This  message  indicates  that  the 
predicted  probability  of  failure  is  greater  than  unity,  which  is  incorrect.  Check  the  soil 
property  input  data  files  for  accuracy.  Take  particular  care  that  the  values  NUM2,  NUM3 
and  NUM4  are  correctly  entered  in  the  SOILDATA.DAT  file,  and  that  the  values  of  the 
probabilities  entered  in  each  of  the  UNITWEIG.DAT,  FRICTION.DAT  and 
COHESION.DAT  files  are  realistic  and  sum  to  unity. 

5.  error  -  computed  probability  is  less  than  0.0.  This  message  indicates  that  the 
predicted  probability  of  failure  is  less  than  zero,  which  is  incorrect.  Check  the  soil  property 


input  data  files  for  accuracy.  Take  particular  care  that  the  values  NUM2,  NUM3  and 
NUM4  are  correctly  entered  in  the  SOILDATA.DAT  file,  and  that  the  values  of  the 
probabilities  entered  in  each  of  the  UNITWEIG.DAT,  FRICTION.DAT  and 
COHESION.DAT  files  are  realistic  and  sum  to  unity. 

6.  error  -  too  much  lateral  erosion.  This  message  indicates  that  the  value  supplied  for 
the  fluvial  erosion  parameter  results  in  so  much  lateral  erosion  that  a  vertical  bank  or 
undercut  profile  develops.  Such  a  profile  cannot  be  analysed  by  this  model,  but  may  be 
analysed  using  Thome  and  Tovey  (1981).  Alternatively,  try  using  a  reduced  value  of 
dW/dZ. 

7.  Type  1  Convergence  failure  at  yf*  =  .  This  warning  message  indicates  that  the 
failure  plane  angle  iteration  procedure  failed  to  converge  to  a  solution,  for  the  failure 
surface  at  the  indicated  dimensionless  bank  height.  If  this  occurs,  the  program  sets  the 
failure  plane  angle  for  that  failure  surface  to  the  default  value  given  by  equation  (40)  in  the 
main  report,  and  computations  continue:  program  execution,  therefore,  does  not  fail.  If  this 
warning  message  is  displayed,  users  should  however  bear  in  mind  the  use  of  the  default 
value  in  interpreting  their  output  data. 

8.  Type  2  Convergence  failure  at  yf*  =  .  This  warning  message  indicates  that  the 
failure  plane  angle  iteration  procedure  converged  to  an  unrealistic  solution  in  which  the 
failure  plane  angle  was  predicted  to  be  negative,  for  the  failure  surface  at  the  indicated 
dimensionless  bank  height.  If  this  occurs,  the  program  sets  the  failure  plane  angle  for  that 
failure  surface  to  the  default  value  given  by  equation  (40)  in  the  main  report,  and 
computations  continue:  program  execution,  therefore,  does  not  fail.  If  this  warning 
message  is  displayed,  users  should  however  bear  in  mind  the  use  of  the  default  value  in 
interpreting  their  output  data. 

9.  Type  3  Convergence  failure  at  yf*  =  .  This  warning  message  indicates  that  the 
failure  plane  angle  iteration  procedure  converged  to  an  unrealistic  solution  in  which  the 
failure  plane  angle  was  predicted  to  be  greater  then  the  bank  angle,  for  the  failure  surface  at 
the  indicated  dimensionless  bank  height.  If  this  occurs,  the  program  sets  the  failure  plane 
angle  for  that  failure  surface  to  the  default  value  given  by  equation  (40)  in  the  main  report, 
and  computations  continue:  program  execution,  therefore,  does  not  fail.  If  this  warning 
message  is  displayed,  users  should  however  bear  in  mind  the  use  of  the  default  value  in 
interpreting  their  output  data. 


E.  TEST  CASE 


This  check  list  example  is  provided  to  ensure  that  the  user  is  familiar  with  the 
procedures  used  to  run  the  model,  and  to  check  that  the  supplied  program  is  working 
accurately.  First  load  the  example  data  files  provided  together  with  an  executable  version  of 
the  source  code  provided  onto  a  computer.  Then  follow  these  steps: 

1.  Type  DARBY  and  hit  the  RETURN  key  to  initiate  program  execution. 

2.  In  response  to  the  prompts  enter  the  following  bank  geometry  values: 

i)  Flood-plain  elevation  (yfp)  =  4.0m 

ii)  Elevation  of  bank  toe  (yf)  =  L5m 

iii)  Bank  angle  =  65  degrees 

iv)  Elevation  of  base  of  bank  slope  (ys)  =  2.2m 

v)  Elevation  of  base  of  vertical  slope  (yt)  =  4.0m 

vi)  Elevation  of  base  of  tension  crack  (yk)  =  3.8m 

3.  Name  the  output  data  file. 

4.  Do  not  select  the  critical  geometry,  pore  water  pressure  or  confining  pressure  analyses. 

5.  Select  the  non-probabilistic  bank  stability  analysis,  check  the  data  and  type  ok  in 
response  to  the  prompts. 

6.  Once  program  execution  is  complete,  check  the  named  output  data  file.  If  the  user 
correctly  specified  the  data  and  options,  and  the  program  is  working  correctly,  the 
following  output  data  should  be  obtained: 


Factor  of  Safety 

2.77 

Failure  Plane  Angle  (degrees) 

52.39 

Failure  Block  Volume  (m''3/m) 

1.40 

Failure  Block  Width  (m) 

1.68 

Rood-Plain  Widening 

.84 

Height  of  Most  Critical  Failure  Surface 

.050 

F.  EXAMPLE  APPLICATIONS 


It  is  recommended  that  the  user  work  through  the  following  examples  in  order  to 
become  familiar  with  the  potential  applications  of  the  model.  Step-by  step  solutions  to  each 
problem  are  provided  below,  but  users  should  attempt  to  work  through  the  examples  on 
their  own  prior  to  comparing  their  results  with  the  given  solutions.  Each  of  the  following 
examples  are  based  on  analysis  of  the  bank  profile  shown  in  Figure  AIVl.  The 
corresponding  bank  material  properties  are  summarized  in  Table  AIVS.  In  the  following 
examples,  the  user  will  apply  the  computer  program  to: 

1.  Determine  the  probability  of  failure  of  the  bank  described  by  Figure  AIVl  and  Table 
AIV5,  without  considering  the  effects  of  pore  water  pressure  or  hydrostatic  confining 
pressure; 

2.  Determine  the  effects  of  including  pore  water  pressure  in  the  analysis  of  the  stability  of 
the  bank,  by  analysing  the  probability  of  failure  for  the  case  of  a  ground  water  table  at 
elevation  4.7  m; 

3.  Determine  the  initial  factor  of  safety  of  the  bank  described  by  Figure  AIV 1  and  Table 
AIV5,  and  the  amount  of  degradation  required  to  erode  the  bank  to  the  point  of  failure, 
assuming  a  fluvial  erosion  coefficient  of  0.1.  For  simplicity,  do  not  analyse  pore  water  or 
hydrostatic  confining  pressure  effects  in  this  example. 

4.  Compare  the  amount  of  degradation  required  to  bring  the  bank  to  failure  when  a  ground 
water  surface  elevation  of  4.7m  is  incorporated  in  the  analysis,  with  the  degradation 
amount  determined  in  example  (3).  Again,  assume  the  fluvial  erosion  coefficient  has  a 
value  of  0.1,  and  exclude  hydrostatic  confining  pressure  effects. 

Example  1:  Solution 

1.  Enter  the  soil  property  data  in  Table  AIV5  to  the  input  files,  UNITWEIG.DAT, 
COHESION.DAT  and  FRICTION.DAT.  Enter  the  modal  values  of  soil  cohesion 
(14,000.0),  friction  angle  (25.0)  and  unit  weight  (23,800.0)  to  the  file  SOILDATA.DAT, 
along  with  the  values  NUM2  (006),  NUM3  (009),  NUM4  (012). 

2.  Use  the  bank  profile  provided  in  Figure  AP/l  to  determine  the  values  of  yfp,  jf,  i,  jg, 
and  Jk.  These  should  be,  respectively,  6.0m,  2.0m,  70.0  degrees,  4.0m,  6.0m  and  5.65m 

3.  Initiate  program  execution  by  typing  DARBY  and  hitting  the  RETURN  key. 
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4.  Follow  the  prompts  to  enter  the  bank  geometry  data  and  name  the  output  data  file. 

5.  Do  not  select  the  critical  geometry,  pore  water  pressure  or  hydrostatic  confining  pressure 
options  in  this  example. 

6.  Select  the  probabilistic  bank  stability  analysis. 

7.  Check  your  input  data  and  re-input  the  data  or  proceed,  as  applicable. 

8.  Wait  for  the  program  to  complete  execution.  If  the  program  successfully  executes,  go  to 
the  next  step.  Otherwise  determine  the  cause  of  the  program  failure,  and  correct  it  by  re¬ 
starting  the  procedure  from  step  1. 

9.  Open  your  named  results  file  and  read  off  the  computed  probability  of  failure.  The 
correct  solution  is  0.645.  If  your  answer  does  not  agree  with  this  value,  check  your  input 
data  carefully. 

Example  2:  Solution 

1.  Enter  the  soil  property  data  in  Table  AIV5  to  the  input  files,  UNITWEIG.DAT, 
COHESION.DAT  and  FRICTION.DAT.  The  correct,  formatted  values  entered  in  each  of 
these  files  are  the  same  as  example  1.  Enter  the  modal  values  of  soil  cohesion  (14,000.0), 
friction  angle  (25.0)  and  unit  weight  (23,800.0)  to  the  file  SOILDATA.DAT,  along  with 
the  values  NUM2  (006),  NUM3  (009),  NUM4  (012). 

2.  Use  the  bank  profile  provided  in  Figure  AIVl  to  determine  the  values  of  yfp,  yf,  i,  y^, 
and  yjj..  These  should  be,  respectively,  6.0m,  2.0m,  70.0  degrees,  4,0m,  6.0m  and  5.65m 

3.  Initiate  program  execution  by  typing  DARBY  and  hitting  the  RETURN  key. 

4.  Follow  the  prompts  to  enter  the  bank  geometry  data  and  name  the  output  data  file  (Give  a 
different  name  from  that  used  for  the  first  example). 

5.  Do  not  select  the  critical  geometry,  or  hydrostatic  confining  pressure  options  in  this 
example.  However,  do  select  the  pore  water  pressure  analysis  and  enter  the  ground  water 
surface  elevation,  GWSE  =  4.7m 

6.  Select  the  probabilistic  bank  stability  analysis. 
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7.  Check  your  input  data  and  re-input  the  data  or  proceed,  as  applicable. 

8.  Wait  for  the  program  to  complete  execution.  If  the  program  successfully  executes,  go  to 
the  next  step.  Otherwise  determine  the  cause  of  the  program  failure,  and  correct  it  by  re¬ 
starting  the  procedure  from  step  1. 

9.  Open  your  named  results  file  and  read  off  the  computed  probability  of  failure.  The 
correct  solution  should  now  be  0.692.  If  your  answer  does  not  agree  with  this  value,  check 
your  input  data  carefully.  Note  that  the  probability  of  failure  computed  in  this  example 
compares  to  the  value  of  computed  in  example  1.  The  effect  of  the  including  the  pore  water 
pressure  arising  from  a  ground  water  surface  elevation  of  4.7m  is,  therefore,  to  increase  the 
probability  of  failure  of  the  bank  from  0.645  to  0.692. 

Example  3:  Solution 

1.  Enter  the  soil  property  data  in  Table  AIV5  to  the  input  files,  UNITWEIG.DAT, 
COHESION.DAT  and  FRICTION.DAT.  The  correct,  formatted  values  entered  in  each  of 
these  files  are  the  same  as  example  1.  Enter  the  modal  values  of  soil  cohesion  (14,000.0), 
friction  angle  (25.0)  and  unit  weight  (23,800.0)  to  the  file  SOILDATA.DAT,  along  with 
the  values  NUM2  (006),  NUM3  (009),  NUM4  (012). 

2.  Use  the  bank  profile  provided  in  Figure  AIVl  to  determine  the  values  of  yfp,  yf,  i,  yg, 
and  y]j.  These  should  be,  respectively,  6.0m,  2.0m,  70.0  degrees,  4.0m,  6.0m  and  5.65m 

3.  Initiate  program  execution  by  typing  DARBY  and  hitting  the  RETURN  key. 

4.  Follow  the  prompts  to  enter  the  bank  geometry  data  and  name  the  output  data  file. 

5.  Select  the  critical  geometry  analysis  and  enter  the  value  of  the  fluvial  erosion  coefficient 
(dW/dZ  =  0.1).  Do  not  select  the  pore  water  pressure  or  hydrostatic  confining  pressure 
options  in  this  example. 

6.  Select  the  non-probabilistic  bank  stability  analysis. 

7.  Check  your  input  data  and  re-input  the  data  or  proceed,  as  applicable. 
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8.  Wait  for  the  program  to  complete  execution.  If  the  program  successfully  executes,  go  to 
the  next  step.  Otherwise  determine  the  cause  of  the  program  failure,  and  correct  it  by  re¬ 
starting  the  procedure  from  step  1. 

9.  Open  your  named  results  file  and  read  off  the  initial  factor  of  safety  and  amount  of 
degradation  required  to  initiate  bank  instability.  The  correct  value  of  initial  factor  of  safety 
is  1.03  -  the  bank  is  predicted  to  be  on  the  verge  of  failure.  The  amount  of  degradation 
required  to  initiate  failure  is  0.10m.  If  your  answers  do  not  agree  with  these  values,  check 
your  input  data  carefully. 

Example  4:  Solution 

1.  Enter  the  soil  property  data  in  Table  AIV5  to  the  input  files,  UNITWEIG.DAT, 
COHESION.DAT  and  FRICTION.DAT.  The  correct,  formatted  values  entered  in  each  of 
these  files  are  the  same  as  example  1.  Enter  the  modal  values  of  soil  cohesion  (14,000.0), 
friction  angle  (25.0)  and  unit  weight  (23,800.0)  to  the  file  SOILDATA.DAT,  along  with 
the  values  NUM2  (006),  NUM3  (009),  NUM4  (012). 

2.  Use  the  bank  profile  provided  in  Figure  AIVl  to  determine  the  values  of  yfp,  yf,  i,  y^,  yj 
and  yjj^.  These  should  be,  respectively,  6.0m,  2.0m,  70.0  degrees,  4.0m,  6.0m  and  5.65m 

3.  Initiate  program  execution  by  typing  DARBY  and  hitting  the  RETURN  key. 

4.  Follow  the  prompts  to  enter  the  bank  geometry  data  and  name  the  output  data  file. 

5.  Select  the  critical  geometry  analysis  and  enter  the  value  of  the  fluvial  erosion  coefficient 
(dW/dZ  -0.1).  This  time,  select  the  pore  water  pressure  analysis  and  enter  the  ground 
water  surface  elevation  value,  GWSE  =  4.7m.  Do  not  select  the  hydrostatic  confining 
pressure  options  in  this  example. 

6.  Select  the  non-probabilistic  bank  stability  analysis. 

7.  Check  your  input  data  and  re-input  the  data  or  proceed,  as  applicable. 

8.  Wait  for  the  program  to  complete  execution.  If  the  program  successfully  executes,  go  to 
the  next  step.  Otherwise  determine  the  cause  of  the  program  failure,  and  correct  it  by  re¬ 
starting  the  procedure  from  step  1. 
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9.  Open  your  named  results  file  and  read  off  the  initial  factor  of  safety  and  amount  of 
degradation  required  to  initiate  bank  instability.  The  correct  value  of  initial  factor  of  safety 
is  0.92.  The  initial  bank  is,  therefore,  predicted  to  be  unstable  with  respect  to  mass  failure 
and  no  further  degradation  is  required  to  initiate  failure.  If  your  answers  do  not  agree  with 
these  values,  check  your  input  data  carefully.  Notice  that  the  effect  of  including  the  pore 
water  pressures  resulting  from  insertion  of  the  ground  water  table  at  4.7m  is  to  decrease  the 
amount  of  degradation  required  to  initiate  failure  from  0.10m  to  Om,  as  compared  to  the 
value  obtained  in  the  third  example.  Notice  also  that  the  initial  factor  of  safety  value  is  also 
reduced,  from  1.03  to  0.92. 

Table  AIV5  Summary  of  soil  property  values  entered  in  innut  data  files  for  examnle 
applications 


UNITWEIG.DAT 

1  COHESION.DAT 

|friction.dat 

P(Y) 

Cohesion 

P(c) 

m 

20,200.0 

0.000 

10,000.0 

0.000 

10.0 

0.100 

21,600.0 

0.150 

11,000.0 

0.000 

15.0 

0.100 

23,800.0 

0.494 

12,000.0 

0.100 

20.0 

0.200 

24,900.0 

0.206 

13,000.0 

0.150 

25.0 

0.200 

26,210.0 

0.150 

14,000.0 

0.220 

30.0 

0.200 

30,000.0 

0.000 

15,000.0 

0.160 

35.0 

0.100 

16,000.0 

0.140 

40.0 

0.100 

17,000.0 

0.120 

45.0 

0.000 

18,000.0 

0.110 

50.0 

0.000 

19,000.0 

0.000 

20,000.0 

0.000 

21,000.0 

0.000 
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